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Preface 

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Calibration and Validation Program has a broad spec- 
trum of responsibilities as outlined in Volume 3 of the SeaWiFS Technical Report Series (McClain et al. 1992a). 
The four primary functions of the program are 1) bio-optical algorithm development, 2) atmospheric correc- 
tion algorithm development, 3) sensor calibration and characterization, and 4) product verification and quality 
control. Each of these categories encompass a number of activities presently underway at the National Aero- 
nautics and Space Administration’s (NASA) Goddard Space Flight Center (GSFC), and within the science 
community as supported by contracts and grants from the SeaWiFS Project and NASA Headquarters. The 
SeaWiFS Project places great emphasis on documentation of its development activities through volumes in the 
SeaWiFS Technical Report Series and the refereed literature, e.g., Hooker et al. 1993. Other volumes of the 
series published under Calibration and Validation Program support include Volumes 4 (McClain et al. 1992b), 
5 (Mueller and Austin 1992), 7 (Darzi 1992), 10 (Woodward et al. 1993), and 14 (Mueller 1993). Because many 
of the studies and other activities undertaken by the Calibration and Validation Program are not extensive 
enough to require dedicated volumes of the series, it was decided to publish volumes composed of brief chapters. 
Volume 13 is the first in a set of such volumes. 

The Calibration and Validation Program relies on the outside research community for the bio-optical and 
atmospheric correction data collection, as well as for algorithm development, but does have the responsibility 
for evaluating and comparing the algorithms and for ensuring that the algorithms are properly implemented 
within the SeaWiFS Data Processing System. Volume 13 consists primarily of sensitivity and algorithm (bio- 
optical, atmospheric correction, and quality control) studies based on the Coastal Zone Color Scanner (CZCS) 
and historical ancillary data undertaken to assist in the development of SeaWiFS specific applications. 

Greenbelt , Maryland — C.R. McClain 

June 1993 
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Abstract 

Although the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Calibration and Validation Program relies on 
the scientific community for the collection of bio-optical and atmospheric correction data as well as for algorithm 
development, it does have the responsibility for evaluating and comparing the algorithms and for ensuring that 
the algorithms are properly implemented within the SeaWiFS Data Processing System. This report consists 
of a series of sensitivity and algorithm (bio-optical, atmospheric correction and quality control) studies based 
on Coastal Zone Color Scanner (CZCS) and historical ancillary data undertaken to assist in the development 
of SeaWiFS specific applications needed for the proper execution of that responsibility. The topics presented 
are as follows: 1) CZCS bio-optical algorithm comparison, 2) SeaWiFS ozone data analysis study, 3) SeaWiFS 
pressure and oxygen absorption study, 4) pixel- by- pixel pressure and ozone correction study for ocean color 
imagery, 5) CZCS overlapping scenes study, 6) a comparison of CZCS and in situ pigment concentrations in the 
Southern Ocean, 7) the generation of ancillary data climatologies, 8) CZCS sensor ringing mask comparison, 
and 9) sun glint flag sensitivity study. 


Prologue 

The purpose of the Sea-viewing Wide Field-of-view Sen- 
sor (SeaWiFS) Project is to obtain valid ocean color data of 
the world ocean for a five-year period, to process that data 
in conjunction with ancillary data to meaningful biologi- 
cal parameters, and to make that data readily available to 
researchers. The National Aeronautics and Space Admin- 
istration’s (NASA) Goddard Space Flight Center (GSFC) 
will develop a data processing and archiving system in 
conjunction with the Earth Observing Satellite Data and 
Information System (EOSDIS), which includes a ground 
receiving system, and will oversee a calibration and vali- 
dation effort to ensure the integrity of the final products. 

The Calibration and Validation Team (CVT) has three 
main tasks: calibration of the SeaWiFS instrument, devel- 
opment and validation of the operational atmospheric cor- 
rection algorithm, and development and validation of the 
derived product algorithms, such as chlorophyll a concen- 
tration. Some of this work will be done internally at GSFC 
while the remainder will be done externally at other insti- 
tutions. NASA and the Project place the highest priority 
on assuring the accuracy of derived water-leaving radiances 
globally, and over the entire mission. If these criteria are 
met, development of global and regional biogeochemical 
algorithms can proceed on many fronts. These activities 
are discussed in detail in The SeaWiFS Calibration and 
Validation Plan (McClain et al. 1992a). 

Because many of the studies and other works under- 
taken with the program are not extensive enough to re- 
quire dedicated volumes of the SeaWiFS Technical Report 
Series , the CVT has decided to publish volumes composed 
of brief, but topically specific, chapters. Volume 13 is the 
first in a set of such volumes and consists primarily of sen- 
sitivity and algorithm (bio-optical, atmospheric correction, 
and quality control) studies based on CZCS and historical 
ancillary data undertaken to assist in the development of 


SeaWiFS specific applications. A short synopsis of each 
chapter is given below. 

L CZCS Bio-optical Algorithm Comparison 

Several pigment concentration algorithms developed for 
use with Coastal Zone Color Scanner (CZCS) data are de- 
scribed and compared. These include the standard two- 
channel switching algorithm of Gordon et al. (1983), an 
iterative algorithm by Smith and Wilson (1981), a three- 
channel algorithm by D. Clark (Muller-Karger et al. 1990) 
and an algorithm developed by the European community 
for the reprocessing of CZCS data from European waters. 

2. SeaWiFS Ozone Data Analysis Study 

Sensitivity analyses are performed on CZCS imagery to 
simulate the impact of erroneous estimates of ozone optical 
thickness on satellite derived water-leaving radiances and 
pigment concentrations. Time series of total ozone from 
the Total Ozone Mapping Spectrometer (TOMS) provide 
an indication of the spatial and temporal variability within 
a CZCS scene. 

3. SeaWiFS Pressure and Oxygen Absorption Study 

Sensitivity analyses are performed on CZCS imagery 
to simulate the impact of erroneous estimates of Rayleigh 
optical thickness on satellite derived water-leaving radi- 
ances and pigment concentrations. Time series of sea level 
pressure from the National Meteorological Center (NMC) 
provide an indication of the spatial find temporal variabil- 
ity within a CZCS scene. 

4. Pixel- by-Pixel Pressure and Ozone Correction Study 

for Ocean Color Imagery 

Sensitivity analyses are performed on CZCS imagery 
to estimate the impact of erroneous estimates of Rayleigh 
and ozone optical thicknesses on satellite derived pigment 
concentrations. 
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5. CZCS Overlapping Scenes Study 

Two coincident scenes from consecutive CZCS orbits 
are analyzed to investigate the consistency of the derived 
products under differing satellite and solar azimuth and 
zenith angles. The Miami DSP CZCS edge mask algorithm 
is described. 

6 . A Comparison of CZCS and In Situ Pigment 

Concentrations in the Southern Ocean 

The large-scale distribution of pigments in the South- 
ern Ocean, as viewed from the CZCS, shows extensive 
blooms and enhanced pigments which are distributed asym- 
metrically about the Antarctic continent. Comparative 
analysis with an extensive database of historical in situ 
data reveals that the magnitude of these enhanced pig- 
ments may actually be 1.8 times higher than previously 
reported. Pigment concentrations are computed using a 
new Southern Ocean CZCS algorithm adjusted to reflect 
regional differences in bio-optical properties of the water 
column and compared to estimates made using an exten- 
sive database of in situ pigment data. 

7. The Generation of Ancillary Data Climatologies 

The method used for generating monthly climatolog- 
ical averages computed on a global basis, of wind speed, 
total ozone, atmospheric pressure, and relative humidity at 


the ocean surface is described. The calibration and valida- 
tion element has computed the climatologies, and placed 
the results in a separate file in the National Center for 
Supercomputing Applications (NCSA) Hierarchical Data 
Format (HDF) for each of the ancillary parameters. This 
chapter describes the method used for the generation of 
the climatologies. 

8. CZCS Sensor Ringing Mask Comparison 

Three different methods for handling CZCS bright tar- 
get recovery are described and compared using a CZCS 
test scene from the Bering Sea. The three techniques are 
the Mueller (Mueller 1988), the SEAPAK (McClain et al. 
1991a and 1991b and Brock et al. 1991) and the Miami 
DSP (Evans and Gordon 1993) methods. 

9. Sun Glint Flag Sensitivity Study 

The statistical wind speed dependent surface slope dis- 
tribution of Cox and Munk (1954a and 1954b) is used to 
estimate the sun glint affected area in a CZCS image. The 
probability of a pixel being contaminated by glitter is a 
function of sea surface wind speed and satellite viewing 
geometry. In the sample case presented, the areal extent 
of the flag expands very little as wind speed increases be- 
yond about 7ms” 1 . 
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Chapter 1 


CZCS Bio-optical Algorithm Comparison 

Charles R. McClain 

NASA Goddard Space Flight Center y Greenbelt, Maryland 
Eueng-nan Yeh 

General Sciences Corporation , Laurel f Maryland 
Abstract 

Several pigment concentration algorithms developed for use with CZCS data are described and compared. 
These include the standard two-channel switching algorithm of Gordon et al. (1983), an iterative algorithm by 
Smith and Wilson (1981), a three-channel algorithm by D. Clark (Muller-Karger et al. 1990) and an algorithm 
developed by the European community for the reprocessing of CZCS data from European waters. The first 
three produce the greatest similarity because they are based on the NIMBUS Experiment Team (NET) data 
set (Clark 1981). 


1.1 INTRODUCTION where m is the index of refraction and p is the Fresnel 

reflectivity. 

The total radiance received by the CZCS is governed 
by the following equation: 1.2 METHODOLOGIES 


Lt( Ai) = t(Xi)Lwpi) + L r (Xi) + L a { Xi) (1) 

where A»=443, 520, 550, and 670 nm, respectively for i=l- 
4; L t is the total radiance; Lw is the water-leaving radi- 
ance; L r is the Rayleigh radiance, L a is the aerosol radi- 
ance, and t is the diffuse transmittance of the atmosphere. 
For a given satellite viewing geometry, L r can be calcu- 
lated from the multiple scattering radiative transfer theory 
(Gordon et al. 1988). L 0 (Aj), is related to L a (670) by the 
expression: 


S( Xi) = 


Lq(Aj) 

La (670) 

*b( Ai)T oz (A i ) 

F 0 (670) T„(670) UV 


( 2 ) 


where n(A*) is conceptually similar to the Angstrom expo- 
nent (Angstrom 1964), Fq is the incident solar irradiance, 
and T oz is the diffuse transmittance of ozone. Making use 
of (2), (1) becomes 


t(Xi)L w (\i) = L t (Xi) - L r ( Xi) - S(Xi)L a (m). (3) 


The water-leaving radiance, Lvv> can be converted to 
subsurface water radiance L s using 


Ls(Xi) 


m 2 Lty( Ai) 
1 ~P 


(4) 


The system described in (3), consists of 4 equations 
(for the first 4 CZCS bands) and 5 unknowns, i.e., Lw(A»), 
i=l-4, and L a (670), if the n(A*) terms are specified. To 
close the system and enable a solution, an additional con- 
dition is required and the various algorithms employ dif- 
ferent approaches in doing so. If the n(Ai) terms are not 
specified, more conditions must be assumed. 

1.2,1 Gordon et al. Method 

Gordon et al. (1983) proposed a total absorption con- 
dition at 670 nm, i.e., L a (670)=0. This is a reasonable as- 
sumption for most low pigment Case 1 (Morel and Prieur 
1977) waters where the optical characteristics are domi- 
nated by phytoplankton and their covarying detrital ma- 
terial. Fig. 1 is a schematic of the Gordon et al. method. 

Usually, the n(A*) terms are assumed to be 0.12, 0.00, 
and 0.00, for 443, 520, and 550 nm, respectively, which are 
typical marine haze values. If the n(Ai) terms are vari- 
able, then assumptions on clear water radiances at 520 
and 550 nm (Gordon and Clark 1981) can be assumed for 
low pigment pixels and the n(Aj) terms can be estimated. 
Examples of studies where the n(Ai) terms were derived 
on a scene-by-scene basis are Barale et al. (1986), Banse 
and McClain (1986), McClain et al. (1988), and Muller- 
Karger et al. (1989). Brock et al. (1991) and Brock and 
McClain (1992) used the same principle to derive summer 
mean n(At) terms for the Arabian Sea. 
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The pigment concentration algorithms utilize ratios of 
L,(Ai) which are interchangable with ratios of L\y(Xi) be- 
cause, the index of refraction and FYesnel reflectivity are 
essentially independent of wavelength over this spectral 
range and cancel out. The bio-optical relation between 
pigment concentration, C, and the subsurface water ra- 
diances of CZCS bands 1 , 2 , and 3 proposed by Gordon 
et al. has two branches where each branch uses an algo- 
rithm based on a ratio of water-leaving radiances, i.e., 
£w(443)/£vv(550) or £w(520)/£w(550). Formulations 
using these bands and radiance ratios are as follows: 


and 


C13 


C23 


1.1298 

3.3266 


Lyy(443) 
Liv (550) 

Lvv(520) 

Lw(SS0) 


-1.71 


-2.40 


(5) 

( 6 ) 


where the C subscripts indicate the bands used. These 
algorithms are developed from field observations from the 
waters of the US East and West Coasts, the Gulf of Mexico, 
and also the Sargasso Sea waters. 

The C 13 algorithm is primarily used for low chlorophyll 
concentration waters (Case 1 ) whereas the C23 algorithm 
encompasses both Case 1 and Case 2 waters. The C 13 
algorithm is used when C 13 is less than 1.5 mg m “ 3 or 
when C 13 is greater than 1.5 mg m -3 , but C 23 is less than 
1.5 mg m~ 3 ; otherwise, C 23 is used. 

In SEAPAK, a PC- based image processing and analysis 
package developed at GSFC (McClain et al. 1991a and 
1991b), if £,(443) is less than 0.15 mW cm -2 /im -1 sr _1 , 
C 23 is used. This is because spurious high pigment values 
can occur when the C 23 switch condition is not satisfied. 
Imperfections in the switching criteria at l.Smgm -3 can 
distort the pigment frequency distribution; a three-channel 
algorithm (Muller-Karger et al. 1990) was developed by D. 
Clark (unpub.) to avoid this type of distortion: 


C = 5.56 


£,(443) + £,(520) 
£,(550) 


n -2.252 


(7) 


1.2.2 Smith and Wilson Method 


Smith and Wilson (1981) proposed an iterative ap- 
proach based on the empirical equation: 


1,(670) = 0.08291£,(443) 


£,( 443 ) ' 

£.( 550 ), 


-1.661 


(8) 


in order to avoid the assumption that £,(670) =0 which 
is particularly unrealistic for coastal waters where pigment 
and sediment concentrations are high (Case 2 waters). In 
the SEAPAK implementation of the algorithm, the pig- 
ment concentrations are calculated using (5) and ( 6 ) rather 
than the pigment algorithms proposed by Smith and Wil- 
son. Their method is illustrated in Fig. 2 . 


1.2.3 European Method 

The so-called European method (Bricaud and Morel 
1987, Andre and Morel 1991, Andersen 1991, and Sturm 
1993) has been developed by the Ocean Colour European 
Archive Network (OCEAN) project, an initiative of the 
Joint Research Center (JRC) of the Commission of the 
European Communities (CEC) and of the European Space 
Agency (ESA), for the purpose of processing CZCS data 
in regions adjacent to the European continent. The at- 
mospheric correction and bio-optical algorithms are sig- 
nificantly different than those described above. Starting 
with ( 1 ), terms are normalized as reflectances, #, and the 
system of equations become: 


Rt( Ai) = Rs(Xi) + Ra(Xi) 

= Rs( At) + G(A i )ft o (670), 


(9) 


where kt = ( R t - Rr)/(qT 2 r ), #a = Ra/{qT 2 T ), and R 8 is 
the subsurface reflectance, T 2 r is the two-way diffuse trans- 
mittance for Rayleigh attenuation, q is the water transmit- 
tance factor, and 


G(Xi) = 


Ra( Ai) 

A a (670) 

T 2r (670) /670V 
T 2r (Ai) \ A / 


( 10 ) 


where 7 is the Angstrom exponent and is independent of 
wavelength. A system of four equations with five unknowns 
is given by (9): #,(A»), t=l-4 and 7 . 

Empirical relations (derived from third-order polyno- 
mial curve fits) between #,(670) and the reflectance ratios, 
#,(443)/#, (550) and #,(520)/#, (550), are introduced to 
complete the nonlinear system. The approach is illustrated 
in Fig. 3. It is iterative, operates on each pixel to deter- 
mine 7 , and, through a sequence of reflectance value tests, 
categorizes the pixel as Case 1 low pigment, Case 1 high 
pigment, or Case 2 . The method assumes the following 
initial conditions: 7 is any specified value, #,(670) = 0, 
and the pixel is assumed to be Case 1 low pigment. The 
iteration incorporates convergence tests for #,(670) and 7 
(Case 1 only). A test for Case 2 water using #,(550) is 
applied and, if satisfied, the calculation branches to the 
Case 2 water algorithm. 

Pigment concentrations are derived from a third order 
polynomial function of #,(443)/#, (550) for Case 1 low 
pigment concentrations. A pigment concentration test for 
low concentration Case 1 water is performed and, if the test 
fails, the calculation branches to the Case 1 high pigment 
concentration algorithm. Once beyond these two tests, an 
Angstrom convergence test is applied. If the convergence 
criterion is satisfied, the final reflectances, pigment con- 
centration, and 7 values are output; if not satisfied, the 
calculation initiates another iteration using the computed 
Angstrom exponent. 
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Fig. 1. Flow diagram for the Gordon et al. (1983) algorithm. 



Fig. 2. Flow diagram for the Smith and Wilson (1981) algorithm. 
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Fig. 4. Histograms among Gordon 2-band, Clark 3-band, Smith-Wilson, and European methods. 


Similar procedures apply to Case 1 high chlorophyll 
pixels except the pigment algorithm which is a function 
of the ratio i? 5 (520)/R s (550). The iteration branches the 
Case 2 algorithm if R 5 (550) exceeds a predefined threshold 
value or the derived pigment is greater than or equal to 
10 mg m -3 . 

For a Case 2 pixel, the subsurface reflectance will be 
estimated using a specified value of the Angstrom expo- 
nent, i.e., no iteration on 7. The iteration scheme is ap- 
plied only using an R s ( 670) convergence test. The pigment 
concentration is derived from either a polynomial relation 
based on i? a (520)/R a (550) or i? a (443)/i? 5 (550) depending 
on whether or not R s ( 443) is greater than the threshold 
value of 0.4. 

1.3 RESULTS 

The algorithms discussed above were implemented in 
SEAPAK and applied to two CZCS scenes. Under these 


comparisons, the Rayleigh optical thicknesses (0.237, 0.123, 
0.098, and 0.044 for 443, 520, 550, and 670 nm, respec- 
tively) were assumed to be constants for both scenes. With 
the exception of the European method, the n(Ai) expo- 
nents were assumed to be constants (0.12, 0.00, and 0.00 
for 443, 520, and 550 nm, respectively). The ozone value 
used was the TOMS data located nearest to the center 
pixel in the scene. These were the standard values used in 
the global CZCS processing (Esaias et al. 1986 and Feld- 
man et al. 1989). 

Plate 1 (orbit 8,889, 28 July 1980, sequential day 210) 
depicts the pigment concentrations derived from the differ- 
ent algorithms. The scene center ozone value is 385 Dob- 
son units (DU). The Angstrom exponent 7 derived from 
the European method is illustrated in Plate 2 which varies 
significantly from location to location. The histograms in 
Fig. 4 and the pigment values in Table 1 indicate the Eu- 
ropean algorithm produces less pigment variability, but 
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Fig. 5. Pigment frequency distributions produced using the Gordon et al. 2-band and European methods for 
scene 2 (Plate 3). 


higher mean pigment values than the other algorithms, 
i.e., approximately 1.75 times higher than the other al- 
gorithms. Also, the Gordon et al. 2-band, Clark 3-band, 
and Smith and Wilson algorithms produced consistently 
similar results. 


Table 1. Pigment mean, p, and standard devia- 
tion, a, for the Norwegian Sea scene (Plate 1, orbit 
number 8,889). The means are in units of mg m“ 3 . 


Value 

Gordon-2 Clark-3 Smith-Wilson European 

<7 

0.717 0.705 0.688 1.248 S 

1.686 1.581 1.191 0.790 


The second scene, Plate 3, is an East Coast scene (orbit 
5,106, 28 October 1979, sequential day 301) with a scene 
center ozone value of 313 DU. Two methods, the Gordon 
2-channel and the European, were compared. As in the 
previous scene, the 7 value (Plate 4) varied significantly 
within the scene. Table 2 compares the standard devia- 
tions and mean pigment concentrations between these two 
methods for the second scene. The European algorithm 
mean pigment value was 29% greater than that derived 
using the Gordon et al. 2-band algorithm while the stan- 
dard deviations between the two were comparable. The 
pigment frequency distributions for the two algorithms are 
shown in Fig. 5. 


Table 2. Pigment mean, /x, and standard devia- 
tion, <7, for the US East Coast scene (Plate 3, orbit 
number 5,106). The means are in units of mg m“ 3 . 


Value 

Gordon- 2 

European 


0.637 

0.823 

<7 

1.822 

1.782 


1.4 CONCLUSIONS 

Four different CZCS pigment algorithms were compared 
in one case and two in a second case. In the former, the 
differences among the Gordon et al. 2-band switching, the 
Clark 3-band non-switching, and the Smith and Wilson 
iterative algorithms produced similar results. This is not 
surprising because all are based on the same bio-optical 
data used for computing pigment concentration. The pig- 
ment concentrations derived from the European method 
averaged about 30-175% higher than the other methods. 
Scene 2 statistics are heavily weighted by the large number 
of low pigment pixels which should produce more favorable 
comparisons because of the relative simplicity of the opti- 
cal properties of the Gulf Stream as compared to scene 1. 
The Angstrom exponents derived using the European al- 
gorithm varied significantly within both scenes with ranges 
greater than would be expected. 
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Chapter 2 


SeaWiFS Ozone Data Analysis Study 

Charles R. McClain 

NASA Goddard Space Flight Center , Greenhelt , Maryland 
Eueng-nan Yeh 

General Sciences Corporation , Laurel f Maryland 
Abstract 

Ozone is a key parameter in the atmospheric correction of ocean color data because visible radiation is dif- 
ferentially absorbed as it passes into and out of the atmosphere. Sensitivity analyses are performed on CZCS 
imagery to simulate the impact of erroneous estimates of ozone optical thickness on the derived water-leaving 
radiances and pigment concentrations. Time series of total ozone from TOMS provide an indication of the spa- 
tial and temporal variability within a CZCS scene. It is concluded that the use of climatologies, e.g., monthly, 
would seriously compromise the primary objective of the SeaWiFS mission, i.e., estimation of surface oceanic 
chlorophyll to within 35%. 


2.1 INTRODUCTION 

The importance of ozone absorption of visible radia- 
tion was recognized early in the development of the CZCS 
atmospheric correction algorithm (Gordon 1978, S0rensen 
1981, Sturm 1981, and Williams et al. 1985). Ozone ab- 
sorbs differentially in the visible and has an absorption 
peak at 602 nm. In the standard CZCS atmospheric cor- 
rection method (Gordon et al. 1983), the total radiance 
received by the CZCS is described by (1). The Rayleigh 
radiance, L r , depends upon the geometric orientation be- 
tween the sun, Earth, and satellite, and, for this analy- 
sis, Lw{ 670) is assumed to be zero. L a (Xi) is related to 
L a (670) through the expression 

La(A<) * (era) * )l “ (670) ’ (11) 

where n(A$) is conceptually similar to the Angstrom expo- 
nent. In the expressions for t , L a (A*) and L r (A*), ozone 
optical thickness appears explicitly. Therefore, errors in 
the estimation of ozone concentration propagate through- 
out all components of the correction in a nonlinear manner. 

As discussed in Williams et al. (1985), S0rensen (1981), 
and Sturm (1981), the original correction scheme assumed 
climatological values of ozone optical thickness based on 
season and latitude range. This procedure was used in 
the early NIMBUS Project processing. However, during 
the global CZCS processing (Esaias et al. 1986 and Feld- 
man et al. 1989), each scene was processed using a single 
ozone concentration value from the gridded TOMS data 


for the same day. The value used was the TOMS value 
located nearest to the center pixel of the first scan line. A 
two-minute CZCS scene covered an area of approximately 
1,600km (along scan) by 800km (along track). While it 
was known that errors in the ozone concentration were 
significant, Andre and Morel (1989) were the first to pub- 
lish sensitivity studies which quantified the magnitude of 
the errors. They found that errors in ozone concentration 
within the natural range of variability (±50 DU) could re- 
sult in errors in the estimated pigment concentration of 
at least 25% at low concentrations and exceed 100% at 
concentrations above 8 mg m -3 . 

The SeaWiFS objectives are to quantify chlorophyll 
a concentrations to within 35% over the range of 0.05- 
50 mg m -3 (Hooker et al. 1992). The work of Andre and 
Morel is strictly theoretical. The purpose of this report is 
to use the operational CZCS algorithm on a sample scene 
to actually show the effects on level-2 products. It must 
be noted that SeaWiFS will use algorithms very similar 
to those developed for the CZCS, but the instrument will 
quantify radiances much more precisely due to its 10-bit 
quantization (CZCS was 8-bit) and higher signal-to-noise 
ratios (greater by approximately a factor of two). 

2.2 DATA ANALYSIS 

The sample scene used was one from the US East Coast: 
orbit 5,106, 28 October 1979, sequential day 301 (McClain 
and Atkinson 1985). This scene includes coastal and open 
ocean water masses providing a wide range of surface pig- 
ment concentrations (approximately 0.1-10 mg m -3 ). The 
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Fig. 6. 1979 time series of daily (solid line) and weekly (dashed line) values of ozone concentration at a point 
(87.5° W,25.5° N) within the CZCS scene. 



Fig. 7. 1979 time series of weekly values of ozone concentration at the four corners of the CZCS scene. 
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RT COMPONENT 

Lt(443) 

Lr(443) 

La(443) 

Lw(443) 

Ll(550) 

Lr(550) 

La(550) 

Lw(550) 

PigxIO 

226 DU 

i m 

6.403 

5.211 

0.376 

1.705 

2.262 

2.089 

0.377 

0.318 

0.640 

301 DU 

wm 

6.403 

5.208 

0.385 

1.693 

2.262 

2.056 

0.380 

0.380 

0.880 

376 DU 

i i 

6.403 

5.205 

0.394 ! 

1.681 

2.262 

2.023 

0.383 

0.442 

1.150 


"nf exponents = (0.12, 0, 0) 
for 443, 520 & 550 nm, respectively 

Fig. 8. Comparison of values of terms in (1) at a low pigment concentration pixel for the three ozone 
concentrations. The data was obtained from orbit 5106 on 28 October 1979. 
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RT COMPONENT 

Lt(520) 

Lr(520) 

La (520) 

Lw(520) 

Lt(550) 

Lr(550) 

La (550) 

Lw(550) 

Pigment 

226 DU mi 

3.207 

2.506 

0.457 

0490 

2.602 

1.918 

0.438 

0 496 

2.362 

301 DU 

3.207 

2.528 

0.451 

0.456 

2.602 

1.950 

0.436 

0.432 

2.915 

376 DU j j 

3.207 

2 550 

0.445 

0.424 

2.602 

1 983 

0.433 

0.369 

3.420 


"n" exponents = (0.12, 0, 0) 
for 443, 520 & 550 nm, respectively 

Fig. 9. Comparison of values of terms in (1) at a high pigment concentration pixel for the three ozone 
concentrations. The data was obtained from orbit 5106 on 28 October 1979. 
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RT COMPONENT 

t(443)* 

t(520>* 

t<56cr 

t(670)* 

LtflO 

u<iO 

La(lc) 

L*hc) 



226 DU 


0.0768 

1.0400 

2.0100 

0.9040 

1.0990 

0.7820 

0.3170 

1.1110 

0.7340 

0.3700 

301 DU 

tm 

0.102 

1.380 

2.680 

1.204 

1.099 

0.776 

0.323 

1.111 

0.740 

0.371 

378 DU 


0.128 

1.730 

3.350 

1.504 

1099 

0.771 

0.328 

1.111 

0.746 

0.365 


* Values x 100 

lc=low concentration, he = high concentration 


Fig. 10. Comparison of ozone optical thicknesses and 670 nm radiance values for the terms in (1). Water- 
leaving radiance is assumed to be zero at 670 nm. The data was obtained from orbit 5106 on 28 October 
1979. 


scene was processed using the atmospheric correction al- 
gorithm of Gordon et al. (1988), the calibration of Evans 
(unpub.) used in the global CZCS processing, and three 
different ozone concentrations (226, 301, and 376 DU). The 
n(Xi) exponents used (0.12, 0.00, and 0.00 for 443, 520, 
find 550 nm, respectively) were those applied to all scenes 
processed during the global CZCS processing. The CZCS 
algorithm requires a fixed set of n(A*) exponents be used 
for all pixels within a scene. Plate 5 provides the pigment 
products from the three analyses and the frequency distri- 
butions of log[pigment]. 

The 301 DU value corresponds to the TOMS value at 
the center of the scene. The other two values are 75% and 
125% of the actual value, respectively. According to Bow- 
man and Krueger (1985), the climatological mean value 
(October 1978 to September 1982) for this site is about 
300 DU with an rms deviation of 30 DU. To examine the 
temporal variability at this location, Fig. 6 provides a one- 
year time series at a location within the scene. As can be 
seen, the daily fluctuations can be as large as 50 DU and 
the annual range is about 100 DU. In order to establish the 
spatial variability within a scene, Fig. 7 provides annual 
time series of weekly ozone concentrations (every seventh 
day) at the four corners of the scene. This information 


indicates the spatial variability can exceed 100 DU. 

To better understand how the terms in (1) vary as 
ozone concentration is varied given constant total radi- 
ances, two locations within the scene are used, one low 
concentration site and one high concentration site. Fig. 8 
presents the low concentration example. Values for 443 nm 
and 550 nm are shown because these two wavelengths are 
used in the pigment concentration algorithm at concentra ir 
tions below 1.5 mg m -3 (Gordon et al. 1983). Note that 
the pigment concentrations are multiplied by 10 in order 
to display the values in the graph. Fig. 9 presents the 
high concentration site values for 520 nm and 550 nm wave- 
lengths because these wavelengths are used for concentra r 
tions above 1.5 mg m“ 3 . Note that the high concentration 
value of approximately 3.0 mg m“ 3 is still low compared to 
the range of values over which SeaWiFS is expected to per- 
form accurately. These examples indicate how the small 
perturbations in the Rayleigh and aerosol radiances pro- 
duce sizable changes in the derived water-leaving radiances 
and pigment concentrations. Fig. 10 provides a compari- 
son of the ozone optical thicknesses and the terms in (1) 
for 670 nm at the low and high concentration pixels. 

According to Bowman and Krueger (1985), the mid- 
latitude locations, such as the present case, have a moder- 
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DAY (10/78 - 9/86) 


Fie. 11. Eight-year time series (October 1978 to September 1986) of ozone concentrations over the Bering 
Sea (170° W, 55° N). 



0 365 730 1095 1460 1825 2190 2555 2920 


DAY (10/78 - 9/86) 

Fig. 12. Eight-year time series (October 1978 to September 1986) of ozone concentrations over the Weddell 
Sea (50° W,65° S). Data gaps are evident in the data and straight lines have been plotted across the gap. 
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ate level of ozone variability with the north polar regions 
having the highest mean and rms deviation values. The 
Bering Sea is one of the world’s most biologically produc- 
tive regions and is of particular interest to US fisheries. 

Fig. 1 1 displays an eight-year time series of ozone con- 
centration over the Bering Sea. The daily fluctuations are 
sizable even when compared to the amplitude of the well 
defined annual cycle. The Southern Ocean will also be of 
great interest because of the plans for Joint Global Ocean 
Flux Study (JGOFS) field studies (Anderson 1992) and 
because of concern about the impact of increased ultravi- 
olet (UV) radiation on phytoplankton in the polar regions 
(Smith et al. 1992). 

Fig. 12 is an eight- year time series of ozone concentra- 
tions over the Weddell Sea. The data from the Weddell 
Sea does not exhibit as strong an annual cycle as does the 
Bering Sea data, but it does indicate a downward trend 
with time. Certainly, studies of the impact of the ozone 
hole on primary production will require accurate estimates 
of zone concentration. 

2.3 CONCLUSIONS 

From the analyses presented, the temporal and spa- 
tial variability of total ozone is sufficiently great to require 
a pixel- by- pixel correction of SeaWiFS data if the chloro- 


phyll accuracy goals of the mission are to be met. Obvi- 
ously, daily global ozone data will not be available at 4 km 
resolution. However, as demonstrated in Fig. 7, the spa- 
tial variability on scales of the order of a swath width is 
great enough to introduce chlorophyll concentration errors 
which exceed the ±35% mission objective. Thus, there 
would be no margin for error due to other error sources 
such as those introduced by deficiencies in the bio-optical 
and aerosol correction algorithms. The present operational 
TOMS grid has variable spatial resolution. It is recom- 
mended that ozone data with the finest spatial resolution 
available, e.g., 5°, be used in the generation of all archived 
SeaWiFS products. To implement a pixel- by-pixel correc- 
tion, an interpolation will be necessary. For quick- look 
processing in near-real time, climatological ozone values 
are perfectly adequate. However, observed values must be 
used for subsequent processings, the first of which will be 
approximately one month after data capture. 

A major concern to the SeaWiFS Project is the avail- 
ability of high quality daily ozone data: the NIMBUS-7 
TOMS is no longer operational, the Pegasus TOMS launch 
has been cancelled, the Meteor TOMS is on an unstable 
platform, and the Upper Atmosphere Research Satellite 
(UARS) ozone products may not be of sufficient quality to 
meet the needs of the Project. 
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Chapter 3 


SeaWiFS Pressure and Oxygen Absorption Study 

Charles R. McClain 
Robert S. Fraser 

NASA Goddard Space Flight Center , Greenbelt, Maryland 
Eueng-nan Yeh 

General Sciences Corporation , Laurel !, Maryland 
Abstract 

Rayleigh scattering is a key factor in the atmospheric correction of ocean color data because visible radiation is 
differentially scattered as it passes into and out of the atmosphere. Sensitivity analyses are performed on CZCS 
imagery to simulate the impact of erroneous estimates of Rayleigh optical thickness on the derived water- leaving 
radiances and pigment concentrations. Time series of sea level pressure from the NMC provide an indication 
of the spatial and temporal variability within a CZCS scene. It is concluded that the use of climatologies, e.g., 
standard atmospheric surface pressure, Pq, would seriously compromise the primary objective of the SeaWiFS 
mission, i.e., estimation of surface oceanic chlorophyll to within 35%. 


3.1 INTRODUCTION 

The removal of the Rayleigh radiance from satellite ob- 
servations of ocean color has always been a critical issue 
due to the strong wavelength dependence of the molecular 
scattering. Early studies (Curran 1972) indicated the fea- 
sibility of estimating ocean pigment concentrations from 
satellite measurements. Initially, the CZCS atmospheric 
correction algorithm was based on a single scattering Ray- 
leigh model and a constant surface pressure, i.e., constant 
Rayleigh optical thicknesses (Gordon 1978, S0rensen 1981, 
Sturm 1981, Gordon et al. 1983, and Williams et al. 1985). 
Subsequent studies investigated and provided correction 
schemes for multiple Rayleigh scattering effects (Gordon 
and Castaho 1987 and Gordon et al. 1988). In the stan- 
dard CZCS atmospheric correction method (Gordon et al. 
1988), the total radiance is given by (1). 

The Rayleigh radiance, which accounts for 80-90% of 
sensor received total radiance, depends upon air molecu- 
lar number density and the orientation between the sun, 
Earth, and satellite and, therefore, is a function of surface 
pressure. Also, for this analysis, Lvv(6 70) is assumed to be 
zero. L a (Xi) is related to L o (670) through the expression 
given in (30). 

In the expressions for t, L a (Aj), and L r (Ai), Rayleigh 
optical thickness, r r , appears explicitly. The optical thick- 
ness is a physical quantity measuring the attenuation power 
of air molecules with respect to a specific wavelength of the 
incident light. The values of r r o, Rayleigh optical thick- 
ness at the standard atmospheric surface pressure, Pq, of 


1,013.25 mb, used in CZCS atmospheric correction compu- 
tations were 0.237 at 443 nm, 0.123 at 520 nm, 0.098 at 
550 nm, 0.044 at 670 nm, and 0.0255 at 750 nm (Gordon et 
al. 1988). At any other surface pressure, 


T r 

L r 


T r 0 



1 - exp (~r r /^) 

1 - exp (-Tro/n) r0 ’ 


( 12 ) 

(13) 


where A P is the pressure deviation and p is the cosine 
of the satellite zenith angle. Errors in the estimation of 
surface pressure propagate throughout all components of 
the correction in a nonlinear manner. 

During the global CZCS processing (Esaias et al. 1986 
and Feldman et al. 1989), each scene was processed using 
a single set of Rayleigh optical thickness values which were 
derived using the standard atmospheric surface pressure, 
P 0 . A two-minute CZCS scene covered an area of approx- 
imately 1,600km (along scan) by 800km (along track). 
While it was known that errors due to variations in the 
surface pressure were significant, Andre and Morel (1989) 
were the first to publish sensitivity studies which quan- 
tified the magnitude of the errors. They found that er- 
rors in atmospheric pressure within the range of variabil- 
ity (±15 mb) could result in errors in the estimated pig- 
ment concentration of at least 10% at low concentrations 
and exceed 100% at concentrations above 10 mg m -3 . The 
SeaWiFS Project objectives are to quantify chlorophyll 
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a concentrations to within 35% over the range of 0.05- 
50 mg m” 3 (Hooker et al. 1992). 

The work of Andre and Morel (1989) is strictly theo- 
retical. The purpose of this report is to use the operational 
CZCS algorithm on a sample scene to actually show the 
effects on level-2 products. It must be noted that SeaWiFS 
will use algorithms very similar to those developed for the 
CZCS, but the instrument will quantify radiances much 
more precisely due to its 10-bit quantitization (CZCS was 
8-bit) and higher signal-to-noise ratios (greater by approx- 
imately a factor of two). 

3.2 DATA ANALYSIS 

The sample scene used was one from the Icelandic wa- 
ters (orbit 9,235, 22 August 1980, sequential day 235), 
which includes coastal and open ocean water masses pro- 
viding a wide range of surface pigment concentrations (ap- 
proximately 0.1-10 mg m“ 3 ). The scene was processed 
using the atmospheric correction algorithm of Gordon et 
al. (1988), the calibration of Evans (unpub.) used in the 
global CZCS processing, and three different surface pres- 
sures (993, 1,013.25, and 1,033.5 mb). The n(A,) exponents 
used (0.12, 0.00, find 0.00 for 443, 520, and 550 nm, respec- 
tively) were those applied to all scenes processed during 
the global CZCS processing. The CZCS algorithm requires 
that a fixed set of Angstrom exponents be used for all pix- 
els within a scene. Plate 6 provides the pigment products 
from the three analyses and the frequency distributions of 
log [pigment]. 

The 1,013.25 mb value corresponds to the standard at- 
mospheric surface pressure used during the global CZCS 
processing. The other two values are 98% and 102% of the 
standard pressure, respectively. To examine the temporal 
variability at this location, Fig. 13 provides a one- year time 
series at a location within the scene. As can be seen, the 
daily fluctuations can be as large as 50 mb and the annual 
range is about 80 mb. 

To better understand how the terms in (1) vary as sur- 
face pressure is varied, given constant total radiances, two 
locations within the scene are used, one low concentration 
site and one high concentration site. Fig. 14 presents the 
low concentration example. Values for 443 and 550 nm 
are shown because these two wavelengths are used in the 
pigment concentration algorithm at concentrations below 
1.5 mg m“ 3 (Gordon et al. 1983). Note that the pigment 
concentrations are multiplied by 10 in order to display the 
values in the graph. Fig. 15 presents the high concentration 
site values for 520 and 550 nm wavelengths because these 
wavelengths are used for concentrations above 1.5 mg m -3 . 
Note that the high concentration value of approximately 
4.0 mg m -3 is still low compared to the range of values 
over which SeaWiFS is expected to perform accurately. 

These examples indicate how small perturbations in 
surface pressure affect the Rayleigh and aerosol radiances 
resulting in sizable changes in the derived water-leaving ra- 
diances and pigment concentrations. Note that in Fig. 15, 


the pigment concentration for the 993 mb case was com- 
puted using the 443 nm water-leaving radiance because of 
algorithm switching. Had 520 nm been selected, the con- 
centration would have been 3.72 mg m” 3 . This points out 
that large errors in the derived pigment can result from a 
combination of surface pressure variability and algorithm 
switching effects. Fig. 16 provides the diffuse transmit- 
tances for the low pigment concentration case and the val- 
ues of the terms in (1) for 670 nm (assuming Lw(670) is 
zero) for both low and high pigment concentration cases. 

3.3 OXYGEN ABSORPTION BAND 

Oxygen in the Earth’s atmosphere moderately absorbs 
radiant energy in a narrow band at 761 nm. Its absorption 
depends on the total optical path length, as determined by 
the solar zenith angle and direction of observation, that the 
solar radiation travels through the atmosphere and on the 
total amount of oxygen, with a 2 a variation of 2.5% in 
amount of oxygen. 

The oxygen absorption band contains 61 absorption 
lines. The transmittance of the band, when computed for a 
narrow band of 6cm” 1 wave numbers or 0.35 nm, is shown 
in Fig. 17. The amount of air that the radiation trav- 
els through is two air masses, which is equivalent to the 
sun being at the zenith, traveling to the surface at a mean 
pressure of 1,013 mb, (shown in the bottom band) and then 
back up through the atmosphere towards the zenith. The 
minimum transmittance is 5% at 760 nm, and a secondary 
minimum of 30% at 763 nm. The band extends from 758- 
770 nm. 

A convenient way of taking the transmittance of the 
oxygen into account for instrument design is to calculate 
the width of an equivalent band that would give complete 
absorption and compare that width with the bandwidth of 
a radiometer. The absorption, A, is given by 

A(k) = 1 -t(k), (14) 

where t(k) is the spectral transmission as a function of 
wave number k (cm -1 ). The equivalent bandwidth, A k, is 
computed by 



where the integration is over the complete wavenumber 
band from &i=12, 970cm” 1 to /c2=13,190cm _1 . Table 3 
gives the equivalent widths for air masses 2 and 3 with a 
surface pressure of 1,013 mb. 


Table 3. Equivalent widths for air masses 2 and 3 
with a surface pressure of 1,013 mb. 


Air Mass 

k [cm x ] 

Width [nm] 

2 

80 

4.6 

3 

89 

5.2 
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Fig. 13. 1980 time series of daily values of NMC surface pressure at a point within the CZCS scene. 
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Lw(443) 
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Lw(550) 
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1013 mb 
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5.262 

4629 

0.197 

0.504 

2.204 

1.816 

0.188 

0.221 

2.770 

1033 mb 

□ 

5.262 

4.709 

0.181 

0.432 

2.204 

1.850 

0.172 

0.201 

3.060 


"n* exponents = (0.12,0,0) 

for 443, 520 and 550 nm, respectively 


Fig. 14. Comparison of values of terms in (1) at a low pigment concentration pixel for the three surface 
pressures. The data was obtained from orbit 9235 on 22 August 1980. 
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Fig. 15. Comparison of values of terms in (1) at a high pigment concentration pixel for the three surface 
pressures. The data was obtained from orbit 9235 on 22 August 1980. 
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Fig. 16. Comparison of the diffuse transmittance values (1) for the low pigment concentration case and 
the values of the 670 nm terms in (1) for both the low and high pigment concentration cases. Water-leaving 
radiance is assumed to be zero at 670 nm. The data was obtained from orbit 9235 on 22 August 1980. 
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Fig. 17 . Oxygen absorption band transmittance for an air mass equal to 2 assuming a resolution of 6 cm 
at 1013.25 mb (bottom curve). 


As an example, consider the radiant energy absorbed 
by oxygen in the SeaWiFS band that extends from 745- 
785 nm. Assuming an air mass of 3, which would occur 
if the zenith angles of the sun and line-of-sight were 48° , 
and a surface pressure of 1,013 mb, Table 3 indicates a 
completely absorbing band 5.2 nm wide would absorb the 
same amount of radiant energy as the entire band. The 
percentage of energy removed from the band, E rem , is 


F*rem — 


5.2 

785 - 745 
13%. 


100 % 


(16) 


Since the limits of surface pressure for SeaWiFS obser- 
vations range from about 985-1,035 mb, the absorption for 
air mass 3 would change only slightly by (1±0.025) 13%. 
If the SeaWiFS zenith angles of sun and line-of-sight in- 
creased from 0° to 60°, the equivalent bandwidth would 
increase from 4.6 to 5.8. 

A precise computation of the oxygen absorption for a 
SeaWiFS band can be done easily with detailed line-by-line 
oxygen absorption spectra. 


3.4 ABSORPTION CORRECTION 


It is found that molecular oxygen has weak absorption 
bands in the red region of the solar spectrum. Comparing 
the equivalent oxygen transmittance bandwidth with the 
radiometer bandwidth, the oxygen optical thickness r ox at 
750 nm can be estimated using concepts from Section 3.3. 
For a plane- parallel, homogeneous absorbing (oxygen at 
750 nm) atmosphere, Beer’s Law states the decrease in the 
radiant intensity is modeled by a simple exponential func- 
tion: 


Ji 


J 0 exp 


Tox 


(17) 


L Mo J 


where Jo is the incident radiant intensity, I\ is the intensity 
after traversing through the absorbing medium, r ox is the 
oxygen optical thickness at 750 nm, and /io Is the cosine of 
the solar zenith angle. 

Similarly, assuming a totally reflecting surface, the re- 
flected radiant energy received by the satellite sensor, I 2 
is _ 


h 


1 1 exp 


“Tox 


(18) 




where fi is the cosine of the satellite zenith angle. The 
absorptivity, defined as the fractional part of the incident 
radiation absorbed by the medium, is given by 


A = 


Jo- 1 2 
Jo ’ 


= 1 — exp 


Mo + M J 


(19) 


The bandwidth of CZCS band 5, 750 nm, is 100 nm. 
Following the equivalent bandwidth concept developed in 
Section 3.3, the absorptivities for 2 and 3 air masses are 


— = 1 — e -2T '“, (20) 

100 v ' 

and 

— = 1 — e _3T,,x , (21) 

100 v ' 

respectively. Solving (20) and (21) for r ox gives an average 

value of 0.02. 

The CZCS atmospheric correction processing uses data 
from the first four bands to derive level-2 products. Band 
5, 750 nm, is used for land and cloud screening only. The 
band 5 radiances were converted to percent albedo, a (Eck- 
stein and Simpson 1991), 


a 


L t { 750) 

t(6)t(0o)F o (75O) 


100 %, 


(22) 
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t(6) = exp 


0.5r r -f r oz -f r ox 
cos 8 


that the number of land or cloud pixels flagged increased 
( 23 ' 1.1% from 165,566 to 167,352. 


where Fo(750) is the incident solar irradiance at 750 nm, 
9 and 0o &re satellite zenith and solar zenith angles, re- 
spectively, and t(9) is the diffuse transmittance (1). All 
pixels with band 5 albedo values exceeding a given thresh- 
old are flagged as land or cloud pixels. Taking the oxygen 
absorption effect into account, the band 5 albedo value will 
increase and more pixels will be identified as land or cloud 
pixels. Results from CZCS Icelandic water scenes show 


3.5 CONCLUSIONS 

From the analyses presented, the temporal and spatial 
variability of surface pressure is sufficiently great to require 
a pixel-by-pixel correction of SeaWiFS data if the chloro- 
phyll accuracy goals of the mission are to be met. Also, 
the oxygen absorption correction is necessary for SeaWiFS 
as the 765 nm band straddles the oxygen absorption band 
and will be used for aerosol corrections. 
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Chapter 4 


Pixel-by-Pixel Pressure and Ozone Correction Study 
for Ocean Color Imagery 

Charles R. McClain 

NASA Goddard Space Flight Center 7 Greenbelt , Maryland 
Eueng-nan Yeh 

General Sciences Corporation Laurel, Maryland 
Abstract 

The global CZCS processing (Feldman et al. 1989) applied constant surface pressure and ozone values to all 
pixels of each two-minute scene. Sensitivity analyses are performed on CZCS imagery to estimate the impact 
of erroneous estimates of Rayleigh and ozone optical thicknesses on the derived pigment concentrations. It is 
concluded that the use of climatologies, i.e., standard atmospheric surface pressure, P 0 > and fixed ozone values 
would seriously compromise the primary objective of the SeaWiFS mission, i.e., estimation of surface oceanic 
chlorophyll a and pigment concentrations to within ±35%. 


4.1 INTRODUCTION 

The retrieval of accurate water-leaving radiances from 
satellite observations requires the computation of ozone 
absorption and Rayleigh scattering effects from two passes 
for both incoming and outgoing light. Previous algorithms 
for the processing of CZCS imagery such as those used 
in the global processing (Esaias et al. 1986; Feldman et 
al. 1989) made several assumptions about the atmosphere. 
One approximation was to use a standard atmospheric sur- 
face pressure, P o =l,013.25 mb, to derive the Rayleigh op- 
tical thickness r r q. These constants (0.237, 0.123, 0.098, 
0.044, and 0.0255 for 443, 520, 550, 670, and 750 nm, re- 
spectively) were applied to all pixels for all CZCS scenes 
regardless of the spatial and temporal variability in sur- 
face pressure. Each scene was also processed using a single 
ozone concentration value which was taken from the grid- 
ded observed data for the same day. 

The ozone value used was the TOMS value located 
nearest to the center pixel in the scene. The n(At) expo- 
nents used (0.12, 0.00, and 0.00 for 443, 520, and 550 nm, 
respectively) were also assumed to be constant for scenes 
and are typical of marine haze aerosols. Note that these 
exponents axe the negative of the n exponents defined in 
Gordon et al. (1983 and 1988). Results from the theoret- 
ical study made by Andre and Morel (1989) indicated er- 
rors due to natural variations in the surface pressure and 
total ozone are significant, i.e., at least 10% at low con- 
centrations to more than 100% at concentrations above 
8 mg m -3 . 


This study extends the pressure and ozone variability 
analyses presented earlier (see Chapters 2 and 3) to include 
pixel- by- pixel corrections. 

4.2 DATA PROCESSING 

Three consecutive scenes in the vicinity of the Green- 
land and Norwegian Seas taken from orbit 9,235 on 22 
August 1980 (or sequential day 235) were used for the 
comparison. The scenes include coastal and open ocean 
water masses providing a wide range of surface pigment 
concentrations (approximately 0.1-10mg m“ 3 ). Plate 7 
illustrates the pigment products from these three (two- 
minute) CZCS scenes (identified on the left of the image) 
as derived using the standard global processing algorithms 
as implemented in SEAPAK (McClain et al. 1991a and 
1991b). The center ozone values for scenes A, B, and C 
were 315, 369, and 361 DU, respectively. 

4.2.1 Pixel-by-Pixel Pressure Corrections 

The surface pressure fields used for this study were de- 
rived from Fleet Numerical Oceanography Center (FNOC) 
gridded data. The FNOC products used are available twice 
per day (0000 and 1200 GMT) and have a 2.5° resolution. 
The data were received from the GSFC Distributed Active 
Archive Center (DAAC), as described by Olsen and Mc- 
Clain (1992), in the NASA common data format (CDF). 
Because of the coarse spatial resolution and the fact that 
the observations are made at local noon, both temporal 
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and spatial interpolations are required. Section 4.3 de- 
scribes the interpolation algorithms. The temporal in- 
terpolation was performed first, then the resultant pres- 
sures were spatially interpolated to the SEAPAK control 
point grid, a 26x26 point navigation array, where the Ray- 
leigh radiance values are computed. The Rayleigh radi- 
ance values at pixels located between the control points are 
subsequently determined by linear interpolation of control 
point Rayleigh radiance values. Using interpolated Ray- 
leigh radiance values is possible because the radiances are 
smoothly varying. 

The control point surface pressure data for scenes A, 
B, and C are listed in Tables 4, 5, and 6, respectively. The 
pressure fields over these three scenes varied from 1,008 mb 
at the right-most edge of scenes A and B to 1,029 mb at 
the bottom left portion of scene A. The overall average of 
1,020 mb is higher than the global CZCS processing value, 
1,013.25 mb. The derived pigment image using the pixel- 
by-pixel pressure corrections is shown in Plate 8. In this 
analysis, a constant ozone value was used for each scene. 
Plate 9 is the difference between the pixel- by- pixel and 
the constant pressure analyses, i.e., Plate 8 minus Plate 7. 
Plate 9 shows the variations in surface pressure within a 
scene can result in a significant change of pigment concen- 
trations. 

4.2.2 Pixel-by-Pixel Ozone Corrections 

The effect of applying a pixel- by- pixel ozone correc- 
tion while assuming a constant pressure was also exam- 
ined. The control point ozone values for scenes A, B, and 
C are listed in Tables 7, 8, and 9, respectively. The ozone 
distributions are quite different from scene to scene. The 
ozone range of scene A was 305-388 DU with values lower 
than the scene center value (315 DU) being located in the 
lower left quadrant of the scene. For scene B, which has 
an ozone range of 302-392 DU, values lower than the scene 
center value (369 DU) dominate the right half of the scene. 
For scene C, which has a range of 321-382 DU, the lower 
left quadrant is occupied by relatively high values (approx- 
imately 361 DU). 

The pigment concentrations resulting from the pixel- 
by-pixel ozone corrections are shown in Plate 10. Plate 11 
presents the pigment difference image computed by sub- 
tracting the constant ozone analysis from the pixel-by-pixel 
analysis (Plate 10 minus Plate 7). Effects of the variability 
of ozone values within a scene are clearly evident. 

4.2.3 Pressure and Ozone Corrections 

The pigment values obtained using the combination of 
pixel- by- pixel pressure and ozone corrections are shown in 
Plate 12. Pigment differences resulting from these correc- 
tions and the global CZCS processing corrections (Plate 12 
minus Plate 7) are shown in Plate 13. Examination of the 
difference image indicates that at 11.14° W, 61.5° N (scene 


A, south of Iceland), for instance, the pigment value in- 
creased 39% from 1.66 mg m~ 3 (Plate 7) to 2.31 mg m~ 3 
(Plate 12). In other locations, the pigment value decreased, 
e.g., 28% from 0.614 mg m~ 3 (Plate 7) to 0.441 mg m“ 3 
(Plate 12) at 10.32° W,70.35° N (scene B, east of Iceland). 
The pigment concentration frequency distributions corre- 
sponding to Plates 7 and 12 are shown in Fig. 18 and the 
frequency distribution of the pigment difference field in 
Plate 13 is shown in Fig. 19. 

4.3 DATA INTERPOLATION 

Values for a desired location and time rarely exist from 
observed or simulated data fields. Therefore, an interpola- 
tion algorithm can be used to estimate the quantities re- 
quired from surrounding values. In the present case, both 
the pressure and ozone are available at predefined time and 
space intervals. 

Assume there are N points surrounding a desired lo- 
cation, and time t is a desired time which is between two 
observation times, t\ and £ 2 - The interpolation scheme 
uses an inverse distance and time weighting algorithm. For 
distance, 


and 


where 


JV 

i= 1 
N 

M« 2 )] = £>K(< 2 ), 

i=1 




1 

dl 


E 


i -1 


dj 


(24) 


(25) 


(26) 


d is the distance from the observation point to the point 
of interest, i is the index of the surrounding point, and by 
definition — 1. 

In the case of time, 


V\ 

t<i-t 

*2-0’ 

(27) 


t - ti 

(28) 

v 2 

<2 — ti 

({vm 

= i/iMti)] + **[V(*2)J, 

(29) 


where the v terms are the temporal weighting factors and 
v\ 4- ^2 — 1- Substituting (24) and (25) into (29), 


(MOD = E^O^i) + ^Vi{t 2 )) 

(30) 

= mm- 

This equation states that the order in which the temporal 
and spatial interpolations are performed does not affect 
the result. 
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Table 4. Deviation of control point surface pressure values from 1,013.25 mb for scene A. Add 1,013.25 mb to the table 

value to compute the actual pressure. 

10 10 10 11 12 12 13 13 13 13 13 13 12 12 12 12 12 11 11 11 9 8 6 4 1 -3 

10 11 11 11 12 12 13 13 13 13 13 13 12 12 12 12 12 12 11 11 10 8 6 4 1 -3 

10 11 11 11 12 12 13 13 13 13 13 13 13 12 12 12 12 12 11 11 10 8 7 4 1 -3 

10 11 12 12 12 13 13 13 13 13 13 13 13 13 12 12 12 12 12 11 10 8 7 5 1 -4 

11 11 12 12 12 13 13 13 13 13 13 13 13 13 12 12 12 12 12 12 11 9 8 5 1 -4 

11 11 12 12 12 13 13 13 13 13 13 13 13 13 12 12 12 12 12 12 11 9 6 5 2 -5 

12 11 12 13 12 13 13 13 13 13 13 13 13 13 13 13 12 13 12 11 11 9 6 6 2 -5 

12 12 12 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 12 11 9 6 6 2 -5 

12 13 12 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 12 12 10 7 6 3 -5 

13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 12 12 10 7 4 3 -5 

13 13 14 13 13 13 13 13 13 13 13 13 13 13 13 13 13 14 14 13 12 10 8 4 3 -4 

13 13 14 14 14 13 13 13 13 13 13 13 13 13 13 13 13 14 14 13 13 11 8 5 3 -4 

13 14 14 14 14 14 14 13 14 13 13 13 13 13 14 14 14 14 14 13 13 11 8 5 2 -3 

13 14 14 14 14 14 14 14 14 14 14 13 13 13 14 14 14 14 14 13 11 11 8 5 1 -3 

14 14 15 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14 13 13 11 11 8 5 1 -3 

14 14 15 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14 13 12 11 8 5 1 -4 

14 15 15 15 15 14 14 14 14 14 14 14 14 14 14 14 14 14 14 13 12 11 9 5 1 -4 

15 15 15 15 15 14 14 14 14 14 14 14 14 14 14 14 14 14 14 13 12 11 9 5 1 -3 

15 15 15 15 15 15 14 14 14 14 14 14 14 14 14 14 14 14 14 13 12 11 9 5 1 -3 

15 15 15 15 15 15 15 14 15 14 14 14 14 14 14 14 14 14 14 14 12 11 9 6 2 -3 

15 15 16 15 15 15 15 15 15 15 15 15 15 15 15 15 15 14 14 14 12 11 9 6 2 -3 

15 16 16 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 14 14 12 10 9 6 2 -3 

15 16 16 16 15 15 15 15 15 15 15 15 15 15 15 15 15 15 14 14 12 10 9 6 3 -2 

15 16 16 16 16 16 15 15 15 15 15 15 15 15 15 15 15 14 14 14 12 10 9 6 3 -2 

15 16 16 16 16 16 16 16 16 15 16 16 15 15 15 15 15 14 14 13 12 10 10 6 3 -2 

15 15 16 16 16 16 16 16 16 16 16 16 15 16 15 15 15 14 14 13 12 10 10 6 3 -2 

Table 5. Deviation of control point surface pressure values from 1,013.25 mb for scene B. Add 1,013.25 mb to the table 

value to compute the actual pressure. 

55555556677788999999887530 
5555555667788899999988753 -1 

5555566677888899999988653 -1 

6555566677888899999987653 -1 

655566777788999 10 10 99987653 -2 

655566777889999 10 10 10 9987651 -2 

666666788889999 10 10 10 9987650 -2 

6666787889899 10 10 10 10 10 9987640 -2 

6 6 6 7 7 8 8 8 8 9 9 9 10 10 10 10 10 10 9 9 8 7 6 4 0 -3 

6 6 6 7 7 8 9 9 9 9 9 10 10 10 10 10 10 9 9 9 7 7 6 3 0 -3 

6 6 7 7 8 8 9 9 10 9 10 10 10 10 10 10 10 10 9 9 7 6 6 3 0 -4 

7 7 7 8 9 9 9 9 10 10 10 10 10 10 10 10 10 10 9 9 7 6 4 2 0 -4 

7 7 7 8 9 10 10 10 10 10 10 11 11 11 11 10 10 10 10 9 7 6 4 2 0 -3 

7 7 8 8 9 10 10 11 10 11 11 11 11 11 11 11 10 10 9 9 8 6 4 2 0 -3 

7 8 8 9 9 11 11 11 11 11 11 11 11 11 11 11 10 10 9 9 8 6 4 2 0 -3 

7 8 8 10 10 11 11 11 11 11 11 11 11 11 11 11 10 10 9 8 7 6 4 2 0 -3 

8 8 8 10 10 11 11 11 12 11 12 11 12 11 11 11 10 10 9 9 7 6 5 2 0 -4 

8 9 9 10 11 11 12 12 12 12 12 12 12 12 11 11 10 10 9 8 8 7 5 2 0 -4 

8 9 9 10 11 11 12 12 12 12 12 12 12 12 11 11 11 10 10 9 8 7 5 3 0 -4 

8 9 10 10 11 11 12 12 12 12 12 12 12 12 11 11 11 10 10 9 8 7 5 3 0 -5 

8 9 10 11 11 12 12 12 12 12 12 12 12 12 12 11 11 10 10 9 9 7 5 3 0 -4 

9 9 10 11 12 12 12 12 12 12 12 12 12 12 12 11 11 10 10 9 8 7 6 3 0 -4 

9 9 10 11 12 12 12 12 12 12 12 12 12 12 12 12 11 10 10 9 8 8 6 3 0 -3 

9 9 10 11 12 12 12 12 12 13 12 13 12 12 12 12 11 11 10 10 8 8 6 4 0 -3 

10 10 10 11 12 12 13 12 13 13 13 13 13 12 12 12 12 11 11 10 9 8 6 4 1 -3 

10 10 10 11 12 12 13 13 13 13 13 13 12 12 12 12 12 11 11 11 9 8 6 4 1 -3 
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Table 6. Deviation of control point surface pressure values from 1,013.25 mb for scene C. Add 1,013.25 mb to the table 
value to compute the actual pressure. 
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Table 7. Deviation of control point ozone value from the value at the center of scene A. Add 315 DU to the table 
value to compute the actual ozone value. 
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Table 8. Deviation of control point ozone value from the value at the center of scene B. Add 369 DU to the table 
value to compute the actual ozone value. 
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Fig. 18. Pigment concentration frequency distributions obtained from the global CZCS and the pixel- by- pixel 
ozone and pressure correction processing schemes shown in Plate 7 (solid line with a mean of 1.67 mg m -3 ) 
and Plate 12 (x with a mean of 1.88 mg m -3 ), respectively. The x-axis is logarithmic. 



- 0.77 - 0.18 0 0.18 0.77 

Pigment Concentrations (mg/m 3 ) 


Fig. 19. The pigment difference frequency distribution for Plate 13. The x-axis has a logarithmic scaling. 


4.4 CONCLUSIONS 

From the analyses presented, the variability of surface 
pressure and total ozone is large enough to produce a rela- 
tive deviation in concentration that can exceed the ±35% 


accuracy goal of SeaWiFS, leaving no margin for error from 
other sources. Climatologies and fixed ozone values are, 
therefore, inadequate, and it is necessary to have pixel- 
by-pixel corrections for these two effects using the highest 
resolution (both in time and space) products available. 
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Chapter 5 


CZCS Overlapping Scenes Study 

Charles R. McClain 

NASA Goddard Space Flight Center , Greenbelt, Maryland 
Eueng-nan Yeh 

General Sciences Corporation, Laurel, Maryland 
Abstract 

Two coincident scenes from consecutive CZCS orbits are analyzed to investigate the consistency of the derived 
products under differing satellite and solar azimuth and zenith angles. General agreement for pigment retrievals 
to within the 35% SeaWiFS accuracy goal was found. Also, the Miami DSP CZCS edge mask algorithm was 
found to work very well and is described. 


5.1 INTRODUCTION 

At high latitudes, swath overlap will provide multiple 
views of cloud free subscenes from consecutive orbits al- 
lowing for greater numbers of samples in the level-3 bins 
for a given data day. In these situations, the products will 
be derived under widely varying solar and spacecraft az- 
imuth and zenith angles. The present study is designed to 
examine the consistency of the derived products using the 
current CZCS atmospheric correction algorithm (Gordon 
et al. 1988). 

5.2 DATA ANALYSIS 

Two CZCS scenes with large cloud free subscenes cov- 
ering the Denmark Strait area were used. Within both 
swaths, two clear subscenes were identified. The two orbits 
are 9,193 (19 August 1980, sequential day 232, 11:05:53.6- 
11:07:53.3) and 9,194 (12:49:23.6-12:51:23.6). The first 
subscene was centered at 17.1° W, 71.0° N near Scoresby- 
sund, Greenland, and the second subscene is centered at 
35.3° W, 65.6° N near Kap Dan, Greenland. From Plate 14, 
the first subscene is near nadir in the 9,193 swath and on 
the right limb of the 9,194 swath. The second subscene is 
on the left limb of the 9,193 swath and near nadir of the 
9,194 swath. Thus, comparisons of data from both limbs 
of the scan can be compared to retrievals near nadir. In 
order to process the scenes using SEAPAK (McClain et 
al. 1991a and 1991b) the original images were subsampled 
using subsampling factors of 3 and 2 for pixel and line 
directions, respectively, to form 512x512 image files. 

The two scenes were processed using the atmospheric 
correction and bio-optical algorithms of Gordon et al. (1983 
and 1988) and the calibration of Evans and Gordon (1993) 


used in the global CZCS processing. The n(Aj) exponents 
used (0.12, 0.00, and 0.00 for 443, 520, and 550 nm, respec- 
tively) were those that were applied to all scenes processed 
during the global CZCS processing. Note that these expo- 
nents are the negative of the n exponents defined in Gordon 
et al. (1983 and 1988). Ozone values derived from the scene 
center location were applied to all pixels within a scene. 
The ozone values used for orbits 9,193 and 9,194 were 350 
and 368 DU, respectively. This difference would not pro- 
duce significant differences in the pigment retrievals. 

In the pigment products shown in Plate 14, a sensor 
ringing mask has been applied using the method devel- 
oped by R. Evans for the CZCS global processing. Also, 
the right-most 41 pixels of each scan line of orbit 9,194 
(hatched area of Plate 14, upper right panel) were also 
excluded from further analysis as this area failed the edge 
mask flag, also developed by R. Evans for the global CZCS 
processing (see Section 5.3). The pigment images were 
remapped to a common universal transverse mercator 
(UTM) projection so pixel-to-pixel comparisons between 
two orbits could be made. Finally, the lower two panels of 
Plate 14 show the remapped images with the ringing and 
edge masks applied. 

Fig. 20 shows the chlorophyll scatterplot (top) and the 
frequency distributions of log[pigment] (bottom) for the 
two scenes. Data values used for Fig. 20 are from only 
unmasked portions of the common areas of the two scenes. 
The bimodal chlorophyll distributions are a result of the 
switching mechanism in the bio-optical algorithm (Muller- 
Karger et al. 1990). 

In the standard CZCS atmospheric correction method, 
the total radiance received by the CZCS is described by (1). 
Nine-pixel average values of terms in (1) and the derived 
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Fig. 20. Chlorophyll scat terp lot (top) and frequency distributions of log [pigment] (bottom) for unmasked 
common pixels of the two pigment scenes. For the latter, orbit 9,193 (solid line) has a mean pigment concen- 
tration of 2.78 mg m~ 3 and orbit 9,194 (+) has a mean pigment concentration of 2.13 mg m -3 . 
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pigment concentrations are listed in Table 10. A summary 
of the spacecraft and solar azimuth and zenith angles is 
provided in Table 11. 


Table 10. Nine-pixel average values of the terms 
in (1) and pigment concentration. 


Parameter 

Value 


Orbit 

9,193 

9,194 

9,193 

9,194 

Latitude [° N] 

71.0 

71.0 

65.6 

65.7 

Longitude [° W] 

17.1 

17.1 

35.3 

35.3 

Term 

Band 

Average Values 

(Lt) 

1 

4.875 

6.192 

6.727 

5.448 


2 

2.691 

3.432 

3.725 

2.967 


3 

2.102 

2.682 

2.888 

2.319 


4 

0.854 

1.113 

1.232 

0.928 

(Lr) 

1 

4.605 

5.732 

6.313 

4.945 


2 

2.418 

2.967 

3.339 

2.556 


3 

1.854 

2.251 

2.543 

1.956 


4 

0.731 

0.888 

1.019 

0.761 

{ La) 

1 

0.149 

0.274 

0.260 

0.201 


2 

0.149 

0.273 

0.258 

0.202 


3 

0.142 

0.258 

0.243 

0.193 


4 

0.123 

0.225 

0.213 

0.167 

(Lw) 

1 

0.137 

0.226 

0.189 

0.344 


2 

0.135 

0.218 

0.146 

0.227 


3 

0.115 

0.198 

0.116 

0.185 


4 

0.000 

0.000 

0.000 

0.000 

(t) 

1 

0.878 

0.822 

0.821 

0.878 

2 

0.919 

0.880 

0.880 

0.918 


3 

0.917 

0.875 

0.876 

0.915 


4 

0.962 

0.942 

0.942 

0.961 

(Chl.) 

0.895 

0.940 

0.519 

0.397 


Table 11. Solar and spacecraft azimuth and zenith 
angles and scan pixel numbers (a full swath has 
1,968 pixels) for the analysis areas of Table 10. 


Parameter 

Value 

Orbit 

9,193 

9,194 

9,193 

9,194 

Latitude [° N] 

71.0 

71.0 

65.6 

65.7 

Longitude [° W] 

17.1 

17.1 

35.3 

35.3 

Pixel Number 

976 

1,818 

114 

966 

Angle 

Angular Values [ 

°] 

Solar Zenith 

61.50 

58.56 

62.63 

55.60 

Satellite Zenith 

23.17 

52.32 

52.75 

23.23 

Solar Azimuth 

125.15 

96.16 

146.58 

118.46 

Satellite Azimuth 

297.72 

207.74 

19.43 

287.84 


5.3 MIAMI EDGE MASK 

In certain situations, the CZCS derived products are 
questionable along the edges of a scan line. R. Evans eval- 
uated the circumstances where pigment retrievals in par- 
ticular are unreasonably high. He found that the poor 
retrievals were a function of tilt angle and developed an 
algorithm to compute a mask for identifying the affected 
area within a scene. Separate algorithms are used for the 
left and right limbs of the scan. 

The EDGEMASK program generates a mask for both sides 
of the scan. For the left limb, a positive integer defines the 
ending pixel number of the left limb mask (the starting 
pixel is always pixel number one) and must be supplied 
by the user. A zero value for the left limb mask indicates 
no mask is applied. Similarly, a positive value can define 
the starting pixel of the right limb. However, when a zero 
value is input for the right limb mask starting pixel, an 
automated algorithm will be used which is based on the 
tilt angle, a. Each CZCS scan line contains 1,968 pixels, 
but the global CZCS processing subsampled the data by 
a factor of four resulting in only 492 pixels. The program 
computes the right edge starting pixel (the ending pixel is 
492) as follows: 

492, 

492 - [3.75(a - 14°)] 

492 - [7.5(a — 18°) 4- 15] , 

462, 

Note, no mask is required for negative tilts. The tilt angles 
for the 9,193 and 9,194 scenes were both +20°. 

5.4 CONCLUSIONS 

Despite large differences in the atmospheric path length 
of the subscenes compared, the pigment products were 
reasonably similar, thus, the atmospheric correction algo- 
rithms are sufficiently accurate to allow the accumulation 
of data into level-3 bins from different orbits. This result 
indicates no selection criterion is necessary for determining 
the best orbit for a particular level-3 bin. Also, the Miami 
edge mask algorithm accurately determined the width of 
the area of questionably high pigment values on the right 
limb of the 9,194 scene. 


if a < 14°; 
if 14° < a < 18°; 
if 18° < a < 20°; 
if 20° < a. 


(31) 
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Chapter 6 


A Comparison of CZCS and In Situ Pigment Concentrations 

in the Southern Ocean 

Kevin R. Arrigo 
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Charles R. McClain 

NASA Goddard Space Flight Center , Greenbelt , Maryland 
James K. Firestone 

General Sciences Corporation, Laurel, Maryland 
Cornelius W. Sullivan 

University of Southern California , Los Angeles , California 
JOSEFINO C. COMISO 

NASA Goddard Space Flight Center , Greenbelt, Maryland 
Abstract 

The large-scale distribution of pigments in the Southern Ocean, as viewed from the CZCS, shows extensive 
blooms and enhanced pigments which are distributed asymmetrically about the Antarctic continent. Compara- 
tive analysis with an extensive database of historical in situ data reveals that the magnitude of these enhanced 
pigments may actually be 1.8 times higher them previously reported. Pigment concentrations computed using 
a new Southern Ocean CZCS algorithm adjusted to reflect regional differences in bio-optical properties of the 
water column agree to within 5% of estimates made using an extensive database of in situ pigment data. This 
result is encouraging and indicates that the unique features of the large-scale data, including the asymmetrical 
distribution around the continent, are real. 


6.1 INTRODUCTION 

Knowledge of global-scale features of phytoplankton 
biomass and productivity is crucial to obtaining a more 
complete understanding of the role of the Southern Ocean 
in the contemporary global carbon cycle of matter (Fasham 
et al. 1990) and the relationship between the distributions 
of primary producers and higher trophic level consumers 
at the basin scale. As the major agents in the sea re- 
sponsible for the transformat ion of approximately 50 x 10 15 
grams of carbon dioxide into fixed organic carbon (in the 
form of phytoplankton biomass), phytoplankton are cen- 
tral to studying and understanding both problems. These 
single celled plants represent a major potential sink for 
atmospheric carbon dioxide in the sea (Platt and Sathyen- 
dranath 1988) as well as the major food of planktonic graz- 
ers (Frost 1991). 

Accurate assessments of the large-scale distribution, 
abundance, productivity, and sedimentation rates of phy- 
toplankton are difficult to obtain by conventional ship- 
based studies with their characteristic low spatial resolu- 


tion. This is especially true for the Southern Ocean where 
severe weather and seasonal coverage by sea ice greatly re- 
stricts access by ships (Sullivan et al. 1988, Comiso et al. 
1990, and Comiso et al. 1993). However, CZCS pigment 
data has been shown to be useful for characterizing surface 
features of the Southern Ocean (McClain et al. 1991c and 
Comiso et al. 1993). 

Satellite data are a very important source of large-scale 
pigment distributions and standing crops at high latitudes, 
but adequate validation is necessary because of the uncer- 
tainties associated with low sun angle, multiple scattering 
effects, and unknown bio-optical characteristics of surface 
waters. Previous quantitative comparisons of in situ pig- 
ments (chlorophyll a plus phaeopigment) with CZCS data 
have been encouraging. Unfortunately, these studies either 
focused on the Northern Hemisphere (Balch et al. 1992) or 
utilized small Southern Ocean data sets which were spa- 
tially restricted (Sullivan et al. 1988, Comiso et al. 1990, 
and Comiso et al. 1993). 

This study is concerned with comparing mean summer 
in situ pigment values south of 30° S latitude with CZCS 
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climatologies averaged over identical temporal and spatial 
scales. The objective of the study is to assess the appli- 
cability of existing pigment algorithms to Southern Ocean 
waters. 


6.2 CZCS METHODS 


An austral summer Southern Ocean (south of 30° S) 
CZCS image was created by averaging monthly (October- 
March) pigment climatologies from 1978 to 1986. All im- 
ages were processed using the standard NASA global pro- 
cessing (GP) algorithm (Gordon et al. 1983) in conjunction 
with the SEAPAK image processing package (McClain et 
al. 1991a and 1991b). Pigment concentrations in excess 
of 10 mg m“ 3 (less than 1% of all samples) were excluded 
from further analyses due to suspected unreliability. Out- 
put of a regional Southern Ocean (SO) algorithm (Mitchell 
and Holm-Hansen 1991) was also approximated by apply- 
ing the appropriate conversion factor to specific pigment 
intervals from the GP algorithm. 

The basic equations for both the GP and the SO algo- 
rithms are: 


log [C + P] 
iog[C 4- P] 


a + b log 


a + b log 


L u (441) 
L u (550) j 

’ L u (520y 
L u (550) J ’ 


for C < 1.5 (32) 
for C > 1.5 (33) 


where [C+P] is mg chlorophyll a plus phaeopigments m -3 
(the former being C and the latter P), and the coefficients 
a and b differ for the GP and SO algorithms as shown 
in Table 12. Note that there is no statistical difference 
between the GP and the SO L u (520)/L u (550) algorithms. 


Table 12, Coefficients for the GP and SO CZCS 
pigment algorithms. The coefficients for the GP 
algorithm are from Gordon et al. (1983) and the 
coefficients for the SO algorithm are from Mitchell 
and Holm-Hansen (1991). 


CZCS Band Ratio and 
Algorithm Coefficients 

Algorithm 

GP SO 

L u (441)/L u (550) a 

0.14 

0.53 

b 

-1.55 

-1.63 

L u (520)/L„(550) a 

0.63 

0.48 

b 

-4.72 

-3.32 


Pigment concentrations were averaged for each image 
a) over the entire region south of 30° S, and b) for each 
1° of latitude between 30° S and 65° S using the program 
GRDMEAN from the SEAPAK image processing package (Mc- 
Clain et al. 1991a and 1991b). The factors used to convert 
GP algorithm output to SO algorithm output are given in 
Table 13. The conversion process requires the GP value be 
multiplied by the listed factor. 


Table 13. L u (440)/L u (560) as seen by the CZCS, 
pigment concentrations are in mg chlorophyll a plus 
phaeopigments m -3 , [C + P], of the GP and SO 
algorithms (see Table 12), and the factor used to 
convert GP pigment to SO pigment. 


L u (440)/L u (560) 

so 

GP 

Factor 

1.00 

3.388 

1.380 

2.455 

1.25 

2.355 

0.977 

2.411 

1.50 

1.750 

0.736 

2.376 

1.75 

1.361 

0.580 

2.347 

2.00 

1.095 

0.471 

2.322 

2.25 

0.904 

0.393 

2.301 

2.50 

0.761 

0.334 

2.281 

2.75 

0.651 

0.288 

2.264 

3.00 

0.565 

0.251 

2.248 

3.25 

0.496 

0.222 

2.234 

3.50 

0.440 

0.198 

2.221 

3.75 

0.393 

0.178 

2.208 

4.00 

0.354 

0.161 

2.197 

4.25 

0.320 

0.147 

2.186 

4.50 

0.292 

0.134 

2.176 

4.75 

0.267 

0.123 

2.167 

5.00 

0.246 

0.114 

2.158 

5.25 

0.227 

0.106 

2.150 

5.50 

0.210 

0.098 

2.142 

5.75 

0.196 

0.092 

2.134 

6.00 

0.183 

0.086 

2.127 

6.25 

0.171 

0.081 

2.120 

6.50 

0.160 

0.076 

2.113 

6.75 

0.151 

0.072 

2.107 

7.00 

0.142 

0.068 

2.101 


6.3 IN SITU METHODS 

Surface (less than 10 m depth) in situ pigment (chloro- 
phyll a plus phaeopigments) data for the Southern Ocean 
(6,183 samples) were extracted from a pigment database 
compiled at GSFC. This in situ database is the most com- 
prehensive set of surface pigment that has thus far been 
assembled for the Southern Ocean. The locations of the 
data stations are shown in Fig. 21. Data were transferred 
to a spreadsheet program and sorted by latitude. Chloro- 
phyll a concentrations were averaged 1) over the entire re- 
gion south of 30° S, and 2) for each 1° of latitude between 
30° S and 65° S. Because many stations did not include 
concurrent phaeopigment concentrations, a chlorophyll a 
to phaeopigment ratio (C/P) of 2.57d=1.51 (iV— 1,070) was 
applied to each computed mean to approximate total pig- 
ment concentrations. To maintain consistency with CZCS 
data screening procedures, all in situ pigment concentra- 
tions in excess of 10 mg m~ 3 (10 out of 6,183 samples) were 
excluded from further analyses. 

Relative frequency distributions for CZCS (as estimated 
by both the GP and the SO algorithms) and in situ pig- 
ments were determined using the program HIST in the 
SEAPAK image processing package (McClain et al. 1991a 
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Fig. 21. Stations where in situ pigment samples were collected 
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— In situ (# of samples) 


(a) 



Chlorophyll a + phaeopigment (^g L' 1 ) 



In situ pigment (pig L' 1 ) 


Fig. 22. a) The total number of data points for CZCS and in situ pigments for each 1° of latitude between 
30°S and 65°S. b) Relative frequency distribution of chlorophyll a plus phaeopigment for the in situ data set 
(N=6,173) and for the CZCS data processed using the SO (Mitchell and Holm-Hansen 1991) algorithm and 
the GP (Gordon et al. 1983) algorithm. Mean in situ and CZCS values are also shown, c) CZCS pigment 
(chlorophyll a plus phaeopigments) concentration, estimated using the GP and the SO algorithms, versus in 
situ pigment concentration. In situ and CZCS pigment data points were calculated as circumglobal means 
for each 1° of latitude between 30°S and 65°S. The line labeled 1:1 is the line of perfect agreement between 
the CZCS and in situ data. 
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and 1991b). In situ pigments were regressed against CZCS 
pigment concentrations as estimated by both the GP and 
the SO algorithms. The region south of 65° S was excluded 
from this analysis due to a scarcity of both CZCS and in 
situ pigment data (Fig. 22a). 

6.4 RESULTS AND DISCUSSION 

South of 30° S, in situ pigment concentrations aver- 
aged 0.58 [C 4- P] (V=6, 173). A similar study (Fukuchi 
1980) reported a mean in situ pigment concentration for 
waters between 35° S and 63° S of 0.38 mg m~ 3 for chloro- 
phyll a alone, or correcting for the presence of phaeopig- 
ments (C/P=2.57il.51 in the Southern Ocean, V=l,070), 

0. 53 [C + F], within 10% of this study’s in situ estimate. 
The mean CZCS pigment concentration computed using 
the GP algorithm currently in wide use (Gordon et al. 
1983) was 0.32 [C + P], suggesting CZCS pigment concen- 
trations previously reported for the Southern Ocean are 
underestimated by approximately 45%. This is a remark- 
able result because it was previously suspected by some 
that high pigment values in the polar regions may be as- 
sociated with retrieval errors due to low solar angles and 
multiple scattering effects. 

The large discrepancy may be attributed to uncharac- 
teristically low pigment-specific absorption and low detri- 
tal concentrations of Southern Ocean waters (Mitchell and 
Holm-Hansen 1991). A recently developed SO pigment al- 
gorithm (Mitchell and Holm-Hansen 1991), which reflects 
these water characteristics, in fact shows considerable im- 
provement. The SO algorithm, which utilizes water leaving 
radiances ( L u ) at 441 and 560 nm yields pigment concen- 
trations that are 2. 1-2. 5 times greater (Table 13) than the 
GP algorithm for L u (AAl)/L u (560) ratios greater than 1, 

1. e., GP [ C 4- P] less than 1.5. The mean summertime 
pigment concentration computed using the SO algorithm 
is 0.55 [C 4- P], within 5% and 4%, respectively, of the in 
situ means reported here and by Japanese investigators 
(Fukuchi 1980). 

The relative frequency distribution of chlorophyll a plus 
phaeopigments (Fig. 22b) for this study’s in situ data set 
also exhibits substantially better agreement with satellite 
observations when the CZCS data are adjusted to fit the 
SO algorithm than when it is processed using the GP algo- 
rithm. The log-normal pigment distribution and the pres- 
ence of some high pigment values (the data set includes 


162 samples, or 2.6%, with [C 4- P] greater than 2 re- 
sults in in situ and CZCS derived mean pigment concentra- 
tions that are substantially higher than mode concentra- 
tions. Fukuchi (1980) reported a similar percentage of in 
situ pigment concentrations (2%) greater than 2 mg m -3 . 
These high pigments were found in the southern-most wa- 
ters south of 63° S, generally consistent with our observa- 
tions from CZCS. 

To further compare the applicability of the new SO al- 
gorithm and the earlier GP algorithm to Southern Ocean 
waters, summer CZCS pigment concentrations derived us- 
ing each algorithm were regressed against similarly com- 
puted mean in situ pigment concentrations (Fig. 22c). 
Both the SO and GP algorithms were able to explain 71% 
of the variability in the data (the regression coefficients are 
the same because the analyses were performed on the same 
data set). This is encouraging and indicates ocean color 
information obtained from satellites is useful for charac- 
terizing pigment concentrations at high latitudes. More- 
over, the best fit slope of 0.91 obtained with the SO algo- 
rithm indicates it was a better approximation of the in situ 
pigment data, than was the standard GP algorithm (0.41 
slope), which was developed from a limited number (49) of 
bio-optical stations in waters around the US. The results 
presented here emphasize the need to develop regional al- 
gorithms in the processing of ocean color data, particularly 
because the next ocean color mission, SeaWiFS, is sched- 
uled for launch in 1994. 
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Abstract 

The SeaWiFS Data Processing System (SDPS) requires climatologies, in the form of monthly averages computed 
on a global basis of total ozone, wind speed, atmospheric pressure, and relative humidity at the ocean surface. 
These data are used during the generation of performance assessment, or quick-look , products within the level-2 
processing stream. The calibration and validation element has computed the climatologies, and placed the 
results in a separate file in the NCSA HDF for each of the ancillary parameters. This chapter describes the 
method used for the generation of the climatologies. 


7.1 INTRODUCTION 

The processing of SeaWiFS data from level- 1 (cali- 
brated radiances) to level-2 (derived products) requires 
four ancillary fields: total ozone, surface values of wind 
speed, atmospheric pressure, and relative humidity, all of 
which are incorporated in the atmospheric correction al- 
gorithm (H. Gordon, personal communication). The Sea- 
WiFS refined products are generated with a 3-week de- 
lay, using ancillary fields subjected to quality control (QC) 
procedures and originally collected in near-real time. How- 
ever, the performance assessment product, since it is gen- 
erated within 24 hours of data collection to provide a quick 
look at the ocean color field, will not have the benefit of 
QC ancillary data and will, therefore, use climatological 
fields (Fig. 23). 

The basic requirements for the data sets used in gener- 
ating the climatologies are that they be available globally, 
at the finest spatial resolution possible, for the longest time 
period possible, and ideally, that related parameters (such 
as the meteorological fields) be found from the same data 
source. This last requirement, if met, would eliminate the 
need for creating a hybrid output where intercomparisons 
between parameters would be made more difficult. 

7.2 METEOROLOGICAL DATA SETS 

In order to choose an appropriate data source for the 
meteorological climatologies, the data archive of the Ocean 
Color Group in the Laboratory for Hydrospheric Processes 
at GSFC was searched. This data archive consists of a 
large group of files in the NASA CDF format (NSSDC 
1991), representing various types of model outputs, cli- 
matologies and field experiment results (Firestone et al. 


1990). Most of the CDF files were originally generated 
by the staff of the NASA Climate Data System (NCDS), 
now the GSFC DAAC (Olsen and McClain 1992). Only 
one data source in the archive, the Comprehensive Ocean- 
Atmosphere Data Set (CO ADS), contained gridded pres- 
sure, humidity, and wind data (Woodruff et al. 1987 and 
Slutz et al. 1985). In addition, COADS data were available 
as monthly means spanning a 45-year period, 1946-1990, 
providing a relatively long time period for the averaging. 
To check whether trends existed in the data collected over 
this period, time series plots were generated at a Northern 
and a Southern Hemisphere site. These results are sum- 
marized in Section 9.3. 

COADS is a cooperative effort involving several agen- 
cies or organizations, including the National Oceanic and 
Atmospheric Administration’s (NOAA) National Climatic 
Data Center (NCDC), NOAA’s Environmental Research 
Laboratories (ERL), the Cooperative Institute for Research 
in Environmental Sciences (CIRES), and the National Cen- 
ter for Atmospheric Research (NCAR). The effort seeks to 
provide a historical record of 70 million surface marine 
data reports starting in 1854. The COADS data in the 
NASA ocean color archive have been quality controlled 
and are organized in decadal groupings of related parame- 
ters ( trimmed groups). Each decade’s gridded product (at 
global 2° latitude by 2° longitude resolution, or 90 lati- 
tude by 180 longitude points) is derived by binning edited 
marine observations, such as those available from ships of 
opportunity, and providing averages for each month within 
the decade (NSSDC Master Directory 1993). As shown in 
Table 14, two COADS trimmed groups from each decade, 
beginning with the 1940s, were needed to compute the cli- 
matology for each meteorological parameter. 
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Fig, 23. Schematic for the retrieval of ancillary data used in the SDPS for the generation of level-2 products. 
Note the use of climatological ancillary fields for the generation of the performance assessment product. 


Table 14. CDF file names and surface meteorological parameters used in climatology creation. 


File Name 

4601-4912_CQADS_MSTG_GR0UP3 
5001-5912_COADS_MSTG_GROUP3 
600 1-69 1 2_C0 ADS-MSTG.GR0UP3 
7001-7912-C0ADS-MSTG-GR0UP3 
C0ADS_MSTG2_GR0UP3_800 1-8912 
C0ADS-G3.90 

460 1-49 12_C0ADS_MSTG_GR0UP4 
500 1-59 12_C0ADS_MSTG_GR0UP4 
6001-69 12 -CO ADS _MSTG_GR0UP4 
7001-7912-COADS-MSTG.GROUP4 
C0ADS_MSTG2_GR0UP4_8001 -89 12 
C0ADS-G4-90 


Parameters 
Relative Humidity 
Relative Humidity 
Relative Humidity 
Relative Humidity 
Relative Humidity 
Relative Humidity 
Wind Speed and Pressure 
Wind Speed and Pressure 
Wind Speed and Pressure 
Wind Speed and Pressure 
Wind Speed and Pressure 
Wind Speed and Pressure 


Temporal Coverage 
January 1946 — ► December 1949 
January 1950 — ► December 1959 
January 1960 — ► December 1969 
January 1970 — ► December 1979 
January 1980 — ► December 1989 
January 1990 — ► December 1990 
January 1946 — ► December 1949 
January 1950 — ► December 1959 
January 1960 — ► December 1969 
January 1970— ► December 1979 
January 1980 — ► December 1989 
January 1990— * December 1990 


7.3 COADS TIME SERIES 

To ascertain whether the COADS monthly data exhib- 
ited any major trends over the January 1946 to Decem- 
ber 1990 sampling period, time series plots were generated 
for two widely separated sites: Melbourne, Australia, at 
150°E,38°S (Figs. 24-26) and Cape Hatteras, North Car- 
olina, at 75°W,35°N (Figs. 27-29). Separate plots were 
run for each of the surface parameters, pressure, wind 
speed, and relative humidity. The data were extracted 
from individual files, one file per decade, in the NASA 
CDF, using the program TIMENV in the VAX SEAPAK 


package developed at GSFC (McClain et al. 1991a and 
1991b). TIMENV produced one ASCII file per site and pa- 
rameter, the contents of which were plotted using Golden 
Software Inc.’s Grapher package on an 80486 PC. 

Figs. 24 and 29 (wind speed) are the only ones appear- 
ing to have a substantial trend. In order to assess whether 
this trend was observed at the global scale, a COADS 
time series of monthly wind speed averages was generated 
(Fig. 30). The series was generated by reading COADS 
data stored in NASA CDF, using an Interactive Data Lan- 
guage (IDL) program. IDL is an integrated environment 
from Research Systems, Inc. providing visualization of sci- 
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Fig. 24. Time series of monthly COADS surface wind speed at Cape Hatteras, North Carolina, 35N,75W, 
1946-1991. The regression for the series is indicated by the solid line. 



MONTHS AFTER JAN. 1946 


Fig. 25. Time series of monthly COADS surface pressure at Cape Hatteras, North Carolina, 35N,75W 
1946-1991. The regression for the series is indicated by the solid line. 
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MONTHS AFTER JAN. 1946 

Fig. 28 . Time series of monthly COADS surface pressure at Melbourne, Australia, 38S,150E, 1946-1991. 
The regression for the series is indicated by the solid line. 



Fig. 29 . Time series of monthly COADS surface relative humidity at Melbourne, Australia, 38S,150E, 
1946-1991. The regression for the series is indicated by the solid line. 


39 



Case Studies for SeaWiFS Calibration and Validation, Part 1 



1940 1950 1960 1970 1980 1990 2000 

YEAR 


Fig. 30. Time series of monthly global average wind speeds from COADS, 1946-1990. 


entific and engineering data arrays, by representing them 
as images or graphical displays (Research Systems, Inc., 
1992a and 1992b). This figure shows a noticeable increase 
of 1-1.5 ms -1 during the sampling period. However, con- 
sistent with the findings of Cardone et al. (1990), it is 
likely that the increase is explained by the increasing use 
of anemometers instead of Beaufort scale in the post-1950 
era. Also, Cardone et al. point out that the anemometers 
measure winds at an average height of 19.3 m, rather than 
the 10 m used for shipborne measurements. In the study, 
the trends are largely eliminated when winds corrected to 
20 m are used throughout. Since the wind (and also, pres- 
sure, humidity, and ozone) climatologies are used primar- 
ily for qualitative comparisons with near-real time data, 
as part of the ancillary QC procedure, the presence of a 
trend would have minimal impact on the generation of Sea- 
WiFS data. Nevertheless, further analyses prior to system 
integration are planned, so as to determine the most ap- 
propriate averaging period for climatology generation. The 
COADS data set remains the best choice for the climatol- 
ogy, since it contains all three meteorological parameters 
with a relatively long time record. 

7.4 OZONE DATA SET 

The total ozone climatology was generated from daily 
ASCII files, covering the period 1 November 1978 to 31 


January 1992, stored on 2 CD-ROMs produced by the 
NIMBUS TOMS Ozone Processing Team (OPT) at GSFC 
(Bowman and Krueger 1985). The data in these files were 
derived from the gridded TOMS (GRIDTOMS) orbital 
data set, and are global with a resolution of 1.0° in lat- 
itude by 1.25° in longitude (180 latitude by 288 longitude 
points). Since the original grid cells were equal in area, the 
OPT used an interpolation scheme poleward of 50° latitude 
to create an equal-angle grid of constant 1.25° longitude 
resolution. 

7.5 CLIMATOLOGY GENERATION 

In order to store and verify the data in the climatolo- 
gies, a set of software programs were written (summarized 
in Table 15). All software was run under the UNIX op- 
erating system on Silicon Graphics, Inc. (SGI) Iris work- 
stations. Program WPHCLM read the COADS GR0UP3 and 
GR0UP4 CDF files from disk, after staging from write-once 
read-many (WORM) times platters attached to a Digital 
Equipment Corp. (DEC) MicroVAX II system and con- 
version to a network CDF form readable under UNIX. 
For humidity climatology creation, GR0UP3 CDF files for 
the decades of 1940-1980 inclusive, and for the year 1990, 
were read and averages, standard deviations, and the num- 
ber of observations were generated by month of the year, 
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e.g., January statistics were computed from January 1946, 
1947, . . . ,1990 data, etc. 

Similarly, GR0UP4 CDF files were read to generate wind 
speed and atmospheric pressure statistics. Program 03CLM 
read daily ozone ASCII files furnished by the OPT for 
the period November 1978 — > January 1992. The files were 
read directly from two CD-ROMs attached to an SGI Iris 
workstation. Averages, standard deviations, and number 
of observations were then generated by month of the year 
and placed in an HDF file (Univ. of Illinois 1989 and 1993), 
using the same specification as for the meteorological fields. 

7.6 DISPLAY AND ANIMATION 

Since the meteorological and ozone climatologies will be 
the only ancillary data source available for generation of 
the SeaWiFS performance assessment products, it was es- 
sential that their accuracy be verified. To accomplish this, 
the IDL software package was used. Version 3.0 of IDL sup- 
ports the reading and writing of files stored in the HDF 
format, for display with the IDL tools. An extensive user- 
created library provides additional capabilities, including 
animation of images used to generate movie loops of the 
monthly statistics in the HDF files. 

Two primary IDL procedures were written to facili- 
tate the visualization of the gridded global climatologies 
stored as Scientific Data Sets (SDS) in HDF files, on an 
X-windows terminal: SDS I MAGE and SDSANIMATE (see Ta- 
ble 15). SDSIMAGE displays the contents of an SDS residing 
in an HDF file to an IDL window of any size and position, 
using user-specified gray scaling, zoom factors, and anno- 
tation. A sample call to SDSIMAGE while under the IDL 
interface is as follows: 

IDL> img = SDSANIMATE (mode, hdfnam, sdslab, 
ref no, flag, xzoom, yzoom, window, 
sclmin, sclmax, xsize, ysize, xpos, 
ypos, title) 

where: 

img A two-dimensional array containing global wind, 
pressure, humidity, or ozone data, in a scaled byte 
(if mode is B or IB) or floating point (if mode is D 
or ID) representation. The size of the dimensions 
of img are governed by the data type (COADS 
or TOMS) and the xzoom and yzoom factors de- 
scribed below. 

mode The mode for running a procedure. If D, return 
floating point data but do not display an image; 
if ID, return floating point data and display a 
scaled image; if B, return scaled byte data but 
do not display an image; and if IB, return scaled 
byte data and also display a scaled image. 

hdfnam HDF name containing the SDS of interest, includ- 
ing full directory specification. 


sdslab SDS object annotation (label) pointing to the 
SDS of interest. If the SDS is located by its ref- 
erence number (ref no), sdslab is ignored. The 
label can be found by running the NCSA Collage 
utility, specifying as input the HDF pointed to by 
hdfnam. 

ref no SDS reference number pointing to SDS of inter- 
est. If the SDS is located by its label, refno is 
ignored. The reference number can be found by 
running the NCSA Collage utility, specifying as 
input the HDF pointed to by hdfnam. 

flag An indicator for whether sdslab or refno will 
be used to identify the SDS of interest. Specify 
f lag=0 to use sdslab, or flag=l to use refno. 

xzoom A zoom factor along the line direction. A bi- 
linear interpolation scheme is used to create val- 
ues on a new grid, given the original grid, xzoom, 
and yzoom. A value of 1 will result in no zooming 
being done. 

yzoom A zoom factor along the pixel direction. A bi- 
linear interpolation scheme is used to create val- 
ues on a new grid, given the original grid, xzoom, 
and yzoom. A value of 1 will result in no zooming 
being done. 

window ID of IDL window where data will be displayed 
(ignored for a mode value of B or D). If nega- 
tive, window will be created, otherwise an exist- 
ing window is assumed. 

sclmin A minimum data value for byte scaling. Data at 
sclmin will be set to 0 gray level (black), with all 
values between sclmin and sclmax being linearly 
scaled to cover the gray range 0-255. 

sclmax A maximum data value for byte scaling. Data 
at sclmax will be set to 255 gray level (white), 
with ail values between sclmin and sclmax being 
linearly scaled to cover the gray range 0-255. 

xsize The size of the IDL window in the x (along-line) 
direction; valid only for an input value of window 
less than 0. 

ysize The size of the IDL window in the y (along-pixel) 
direction; valid only for an input value of window 
less than 0. 

xpos The x position of the lower left corner of the IDL 
window created (in device coordinates). 

ypos The y position of the lower left corner of the IDL 
window created (in device coordinates). 

title The title of the IDL display window. 

The second procedure, SDSANIMATE, makes successive 
calls to SDSIMAGE, one per month for the parameter of in- 
terest. SDSANIMATE displays a histogram and scaled image 
for each month in the climatology, at the original grid’s 
resolution. The window with each scaled image exactly 
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Fig. 31. Comparison of single month ozone distributions versus 13-year climatology. 


overlays a window containing a histogram for the image, 
so the user can switch back and forth between the windows. 
After the 12 monthly climatology images and histograms 
are displayed, a loop of the images is begun in a seaparate 
window. The user can interactively change the looping 
speed or the color palette applied to the loop. 

7.7 DISCUSSION 

SDSANIMATE was run on each of the four parameter cli- 
matology HDF files as a means of checking the data. The 
results for the monthly average fields are shown in Plates 
15-18. Plate 15 shows the monthly average ozone based on 
data collected from 1978-1991. Note the migration of the 
black bands near the poles, indicating a lack of data due to 
the low sun angle during the winter season in each hemi- 
sphere. Also note the low ozone values at the South Pole 
in October, consistent with findings described in the liter- 
ature. Likewise, Plates 16, 17, and 18 illustrate familiar 
global meteorological patterns, such as the zones of high 
pressure, low relative humidities, and light winds found 
over the subtropical oceans, as well as, the expected sea- 
sonal fluctuations in the magnitude and location of these 
and other global circulation features. 

During the course of this analysis, it has become clear 
that the averaging period of 13 years for the ozone clima- 
tology is too long, given the downward trend in average 


global ozone during this time. It should also be noted that 
the downward trend for 1992 and 1993 was even more pre- 
cipitous than for the period 1978-1991 (Fig. 31). Fig. 31 
shows the percentages of observations within each month 
which are more than 1 standard deviation (lcr) above the 
13- year climatology standard deviation for that month, as 
compared with the percentages of observations more than 
lcr below the 13-year climatology standard deviation. The 
percentages are plotted for each month during the averag- 
ing period. Note the general tendency for more observa- 
tions to have negative percentages over time, that is, there 
are relatively more observations one or more standard de- 
viations below the climatological value than there are one 
or more standard deviations above climatology. 

The trend for more negative percentages has become 
even more dramatic since 1991, perhaps due to the ef- 
fects on atmospheric chemistry and circulation of the Mt. 
Pinatubo volcanic eruption in the Philippines (Kerr 1993 
and Gleason et al. 1993). A full 40% more of the observa- 
tions are greater than lcr below the mean than are greater 
than lcr above the mean, indicating a highly shifted distri- 
bution. Therefore, it is likely that a new ozone climatology 
will be generated prior to the launch of SeaWiFS for use 
in the SDPS, covering a much shorter period on the order 
of 2-3 years. If ozone continues on a downward trend over 
time, the climatology would need to be computed period- 
ically during the lifetime of SeaWiFS. 
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Chapter 8 


CZCS Sensor Ringing Mask Comparison 

Charles R. McClain 

NASA Goddard Space Flight Center, Greenbelt, Maryland 
Eueng-nan Yeh 

General Sciences Corporation, Laurel , Maryland 
Abstract 

Three different methods for handling CZCS bright target recovery are described and compared. The three 
techniques are the Mueller (Mueller 1988), SEAPAK (McClain et al. 1991a and 1991b, and Brock et al. 1991), 
and Miami DSP (Evans and Gordon 1993) methods. A CZCS test scene from the Bering Sea (orbit 2,746 from 
10 May 1979) which includes both large and broken clouds, is used for the comparison. The three methods yield 
similar results in terms of the number of pixels flagged as bright target contaminated and the average pigment 
values within the scene. If the additional mask value in the Miami DSP method is set to the minimum value 
of one pixel, however, the total number of ringing pixels decreases substantially. The SEAPAK analysis using 
the Brock et al. inputs seriously underestimates the ringing in the test scene. This is not surprising because 
the Brock et al. input parameters are optimized for a desert- to-ocean transition rather than a cloud-to-ocean 
transition. 


8.1 INTRODUCTION 

The CZCS bright target recovery, also referred to as 
ringing and electronic overshoot, occurred on the down- 
scan side of bright objects, such as clouds, ice, or desert 
regions onto darker areas. The result was anomalous data 
that could not be corrected resulting in erroneous derived 
products. The bright target recovery characteristics of the 
CZCS were not well characterized prior to the launch of 
the sensor. The occurrence of CZCS ringing was some- 
what erratic, at times difficult to predict, and the down- 
scan extent of the effect was frequently irregular. As a 
consequence, different analysis methodologies by different 
investigators were developed in an attempt to determine 
where and to what (down-scan) extent the data were cor- 
rupted. 

Three different techniques to quantify CZCS ringing, 
i.e. , the Mueller (1988), Miami DSP, and SEAPAK, (Mc- 
Clain et al. 1991a and 1991b, and Brock et al. 1991) meth- 
ods, are described and compared. This study was under- 
taken in preparation for the launch of SeaWiFS which will 
also exhibit bright target recovery characteristics. It is 
expected, however, that the effect will be better quanti- 
fied as a function of bright target exposure duration and 
radiance during the prelaunch sensor characterization by 
Hughes/Santa Barbara Research Center (SBRC), the com- 
pany responsible for building the SeaWiFS instrument. 


8.2 THE SEAPAK METHOD 

The SEAPAK method utilizes Level-1 670 nm (band 4) 
and 750 nm (band 5) data. Band 5 is used to identify land 
and cloud pixels and band 4 is used determine if a ringing 
mask is to be applied. Fig. 32 shows a schematic of the 
algorithm. Values for band 4 (SATGREY), band 5 (LANCLD) 
and pixel-to-pixel band 4 difference (DELTA) thresholds, 
and the number of pixels (equated to DISTANCE) to be 
masked in the down-scan direction are adjustable on a 
scene- by-scene basis. If a pixel value exceeds the band 4 
and band 5 threshold values, the band 4 grey level value is 
compared to the value of the adjacent down-scan pixel and, 
if the difference exceeds the value of DELTA, the subsequent 
N pixels (where N is equal to DISTANCE) are masked. 

The algorithm is not computationally intensive and is 
implemented in both user interactive and batch processing 
modes. The interactive mode allows the input parameters 
to be varied and the product to be viewed without exit- 
ing the program, thereby allowing fine tuning of the input 
parameter values. 

8.3 THE MUELLER METHOD 

The Mueller (1988) method uses only band 5. Unlike 
the SEAPAK method, the ringing distance (in number of 
pixels) varies as a function of the cumulative excess target 
brightness and the sensor gain factor G (G equals 1.00, 
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Fig. 32. Processing schematic for the SEAPAK method. A pass means the pixel value exceeds the threshold 
value. 


1.25, 1.50, and 2.10 for gains 1-4, respectively). The excess 
target radiance, B, is computed as: 


B = L t { 750) 
B = 0, 


Lc{ 750) if L t (750) > 

G ’ 

if L, (750) < 


L c (750) 

G 

L c ( 750) 

G 


(34) 


where L*(750) is the toted 750 nm radiance and L c (750) is 
the cloud radiance threshold (2.45 mW cm -2 pm -1 sr -1 ). 
The effect of a bright pixel is assumed to be negligible for 
distances greater than 50 pixels. 

The B values are averaged over 10 pixels to generate 
five (B) values which are used to compute DISTANCE using 
the formulation 


DISTANCE = a + /?[ln<B), + lnG], (35) 

where a and (3 were estimated to be 3.9 and 30.8, respec- 
tively, and 

(B) a = X>) t e-°- 32 \ (36) 

1 = 1 

The input values for the Mueller method used on the test 
scene are listed in Table 16. 


Table 16. Algorithm input values used for this 
study. 


SEAPAK Method 

LANCLDf 

Grey Level = 21 

SATGREY 

Grey Level = 210 

DELTA 

Grey Level = 10 

DISTANCE 

4 and 20 pixels 

Mueller Method 

a (intercept) 

3.9 pixels 

(3 (slope) 

30.8 pixels 

L c (750) 

2.45 mW cm -2 pm -1 sr~ 1 

Miami DSP Method 


Saturation Threshold 
Difference Threshold 
Cloud Mask Valued 
Maximum Distance 
Additional Mask 


Grey Level = 60 
Grey Level = 9 
Grey Level = 255 
40 pixels 
20 pixels 


fBand 5 


fCorresponds to LANCLD— 21 


8.4 THE MIAMI DSP METHOD 

The method applied in the Miami DSP system uses 
the CZCS level-2 normalized water-leaving radiances at 
520 nm, Lww(520), as a reference to determine the extent 
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of sensor overshoot in the open ocean (Evans and Gor- 
don 1993). In the generation of the level-2 products, a 
land and cloud mask is applied using a constant thresh- 
old radiance in band 5. Pixels down-scan of a bright tar- 
get are flagged if Lvw(520) exceeds a saturation threshold 
(0.7 mW cm -2 fim~ l sr -1 ) or the difference between suc- 
cessive pixels exceeds a difference threshold (0.1 mW cm -2 
/im -1 sr” 1 ). The normalized water- leaving radiance val- 
ues for the range of 0-3.0 mW cm -2 /rm -1 sr _1 are scaled 
to grey scales from 0-255, respectively. 

The values of 0.7 and 0.1 mW cm -2 p" 1 sr“ 1 equate 
to grey levels of 60 and 9, respectively, for the saturation 
and difference thresholds. The threshold tests are applied 
over a maximum distance of 40 pixels in down-scan direc- 
tion of a bright pixel. Finally, additional pixels are added 
unconditionally to the ringing mask. For the global CZCS 
processing (Feldman et al. 1989), 20 pixels (or equivalently, 
5 GAC pixels) were added to the mask. Input values for 
the Miami DSP method are listed in Table 16 and are the 
equivalent values as those used in the global CZCS pro- 
cessing. 

8.5 RESULTS 

The CZCS scene used for this study is a Bering Sea 
scene (orbit 2,746, 10 May 1979, and sequential day 130) 
which was chosen because it includes both large clouds and 
broken clouds. Plate 19 depicts the pigment concentra- 
tions with white (a grey value of 255) representing pixels 
flagged by the land/cloud and ringing masks. The upper 
left panel represents the case where only the land/cloud 
mask (band 5 at 21 counts) was applied resulting in 45% 
of the data being identified as either land or clouds. Pixels 
seriously affected by ringing can be easily identified as very 
low pigment on the down-scan side of the clouds. The re- 
sults of the three ringing mask algorithms are summarized 
in Tables 17 and 18. 

The SEAPAK input values applied to the Bering Sea 
test scene (Plate 19) are listed in Table 16 for two cases. 
The first case applies values derived by Brock et al. 


(1991) based on analyses of many scenes from the Ara- 
bian Sea (DISTANCE equal to 4) . Their analyses emphasized 
data from scans transitioning from the Arabian desert to 
the ocean. The second case uses DISTANCE equal to 20 
based on an interactive analysis of the Bering Sea scene. 

Table 17. Results of sensor ringing mask algorithm 
comparison. The total number of pixels in the scene 
is 262,144. For SEAPAK 1 DISTANCE is equal to 4 
and for SEAPAK 2 DISTANCE is equal to 20. 


Method 

Used 

Flagged 

(F) 

Ringing 

(R) 

R/^NoMask 

Ratio 

No Mask 

117,123+ 

0 

0.0% 

SEAPAK 1 

144,531 

27,408 

23.4 

SEAPAK 2 

196,145 

79,022 

67.5 

Mueller 

212,829 

95,706 

81.7 

Miami DSP 

210,413 

93,290 

79.7 


f The total number of land and cloud pixels. 


Table 18. Mean (/i) pigment concentrations and 
standard deviations (a) for various ringing mask al- 
gorithms (valid range is 0.042-37.8 mg m~ 3 ). 


Method 


o 

No Mask 

2.40 

4.14 

SEAPAK 2 

3.50 

4.94 

Mueller 

3.36 

4.79 

Miami DSP 

3.74 

5.25 


The SEAPAK (DISTANCE equal to 20), Mueller, and 
Miami DSP methods yielded similar results in terms of 
the number of pixels flagged and average pigment values. 
If the additional mask value in the Miami DSP method is 
set to the minimum value of one pixel, the total number 
of ringing pixels decreases substantially to 51,720. The 
SEAPAK analysis using the Brock et al. inputs seriously 
underestimated the ringing in the Bering Sea scene (not 
shown). This is not entirely surprising because the param- 
eters were optimized for a desert- to-ocean transition rather 
than a cloud-to-ocean transition. 
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Chapter 9 


Sun Glint Flag Sensitivity Study 

Charles R. McClain 

NASA Goddard Space Flight Center , Greenbelt , Maryland 
Eueng-nan Yeh 

General Sciences Corporation, Laurel, Maryland 
Abstract 

The statistical wind speed dependent surface slope distribution of Cox and Munk (1954a and 1954b) is used to 
estimate the sun glint affected area in a CZCS image. The probability of a pixel being contaminated by glitter 
is a function of sea surface wind speed and satellite viewing geometry. In the sample case presented, the areal 
extent of the flag expands very little as wind speed increases beyond about 7 m s“ 1 . 


9.1 INTRODUCTION 

The ocean surface waves deform the ocean surface re* 
suiting in randomly oriented facets which reflect incoming 
solar radiance in many directions other than the direction 
of specular reflection from a flat surface. Directly reflected 
light can result in saturation of the CZCS bands and the 
standard atmospheric correction algorithm (Gordon et al. 
1988) does not explicitly estimate the sun glint radiance, 
although, low levels of sun glint are removed from the im- 
agery because the algorithm interpreted glint radiance as 
aerosol radiance. Ocean color sensors are equipped with 
the capability of tilting the scan plane away from the sun 
in order to minimize the glint effect. However, under cer- 
tain conditions (low solar zenith angle or moderate to high 
wind speeds) it is difficult to avoid sun glint. 

The statistical model of wind wave slope distribution 
proposed by Cox and Munk (1954a and 1954b) and dis- 
cussed in Viollier et al. (1980) assumes that the distribu- 
tion of wave facet orientation is approximately isotropic 
and Gaussian and is solely a function of surface wind ve- 
locity. Consequently, with knowledge of the solar and ob- 
servation viewing geometry, the pattern and intensity of 
the reflected solar radiance can be estimated. 

9.2 METHOD 

Plate 20 depicts pigment concentrations derived from 
a CZCS Mediterranean Sea image (orbit 2,573, 20 April 
1979, sequential day 110) having a scanner tilt angle of 
zero degrees. The relatively high pigment values (blue and 
green areas of upper left panel) near the scene center are 
the result of sun glint corruption, and should be flagged as 


invalid. The CZCS level-2 processing and glint flag analy- 
ses processing were accomplished using SEAPAK (McClain 
et al. 1991a and 1991b). 

The probability of a pixel being contaminated by glitter 
is a function of sea surface wind speed, W , and satellite 
viewing geometry, namely, the solar azimuth and zenith 
angles ($0 and Go, respectively) and satellite azimuth and 
zenith angles ($> and 6 , respectively). A probability pa- 
rameter, Per, is defined by 


Prr = 


7TCT* 


exp 


— tan 2 0 n 


(37) 


where a 2 is the mean square surface slope distribution 
which increases linearly with wind speed: 

a 2 = 0.003 + 0.00512 W (38). 


The zenith angle, 0 n , of the vector normal to the surface 
vector for which glint will be observed can be derived from 
the surface reflection angle, uj: 


cos 2u; 
where 


cos6 cos 60 -f sin 6 sin 9q cos(4> — $ 0)1 (39) 


cos 


cos 9 + cos #o 
2cosu 


(40) 


A pixel will be flagged as sun glint contaminated if the 
calculated value of P a is greater than, or equal to, a given 
threshold value. For a known viewing geometry, the sea 
surface wind speed and the assigned threshold value, (an 
estimate of the sun glint affected area) can be determined. 
The number of pixels flagged decreases as the threshold 
value increases, and increases with wind speed (37). 
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Plate 20 shows the sun glint flags for different sea sur- 
face wind speeds (1, 3, 5, 7, 9, 11, and 13ms" 1 ) and 
threshold values (1.0, 1.5, and 2.0). The First GARP 
(Global Atmospheric Research Program) Global Experi- 
ment (FGGE) 1,000 mb wind speed and the FNOC surface 
wind speed at the center of the scene at local noon were 
both 5ms" 1 . This indicates that a P a value between 1.5 
and 2.0 is optimal, i.e., all contaminated pixels are flagged 
with a minimum amount of loss of valid data. 

9.3 CONCLUSIONS 

If the ocean were a perfectly flat surface, the mirror im- 
age of the reflected sun would be seen at the single specular 


point. As the wind speed increases, the ocean surface is 
ruffled by the wind and the glint area expands. The Cox 
and Munk algorithm can be used to estimate the areal ex- 
tent of the contaminated region. The procedure described 
above requires values for P 0 and W. In the case presented, 
a value of P a = 1.5 would be an adequate value. This ex- 
ample indicates that the areal extent of the flag expands 
very little as wind speed increases beyond about 7ms" 1 . 
Other approaches would be to estimate the sun glint radi- 
ance as described in Viollier et al. (1980) and remove the 
glint explicitly from the total observed radiance, or simply 
to set a flag based on some threshold value of the glint 
radiance. 
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ASCII 


Case 1 
Case 2 


CDF 

CD-ROM 

CEC 

CIRES 

COADS 

CVT 

CZCS 

DA AC 
DEC 
DSP 


EOSDIS 

ERL 

ESA 

FGGE 

FNOC 

GARP 

GMT 

GP 

GRIDTOMS 

GSFC 

HDF 

IDL 

JGOFS 

JRC 

Level-0 
Level- 1 
Level-2 
Level-3 

NASA 

NCAR 

NCDC 

NCDS 

NCSA 

NET 

NIMBUS 


NMC 

NOAA 

NSSDC 

OCEAN 

OPT 

PC 

QC 

RSMAS 


Glossary 

American Standard Code for Information Inter- 
change 

Water whose reflectance is determined solely by ab- 
sorption. 

Water whose reflectance is significantly influenced 
by scattering. 

(NASA) Common Data Format 
Compact Disk-Read Only Memory 
Commission of the European Communities 
Cooperative Institute for Research in Environmen- 
tal Sciences 

Comprehensive Ocean- Atmosphere Data Set 
Calibration and Validation Team 
Coastal Zone Color Scanner 

Distributed Active Archive Center 

Digital Equipment Corporation 

Not an acronym, an image display and analysis 

package developed at RSMAS University of Miami. 

Earth Observing System Data and Information Sys- 
tem 

(NOAA) Environmental Research Laboratories 
European Space Agency 

First GARP Global Experiment 
Fleet Numerical Oceanography Center 

Global Atmospheric Research Program 
Greenwich Mean Time 
Global Processing (algorithm) 

Gridded TOMS (data set) 

Goddard Space Flight Center 

Hierarchical Data Format 

Interactive Data Language 

Joint Globed Ocean Flux Study 
Joint Research Center 

Raw data. 

Calibrated radiances. 

Derived products. 

Gridded and averaged derived products. 

National Aeronautics and Space Administration 
National Center for Atmospheric Research 
(NOAA) National Climatic Data Center 
NASA Climate Data System 
National Center for Supercomputing Applications 
NIMBUS Experiment Team 

Not an acronym, a series of NASA experimental 
weather satellites containing a wide variety of at- 
mosphere, ice, and ocean sensors. 

National Meteorological Center 

National Oceanic and Atmospheric Administration 

National Space Science Data Center 

Ocean Colour European Archive Network 
Ozone Processing Team 

(IBM) Personal Computer 

Quality Control 

Rosenstiel School for Marine and Atmospheric Sci- 
ences (University of Miami) 


SDPS SeaWiFS Data Processing System 
SDS Scientific Data Set 

SEAPAK Not an acronym, an image display and analysis 
package developed at GSFC. 

SeaWiFS Sea-viewing Wide Field-of-view Sensor 
SGI Silicon Graphics, Inc. 

SO Southern Ocean (algorithm) 

TOMS Total Ozone Mapping Spectrometer 

UARS Upper Atmosphere Research Satellite 
UTM Universal Transverse Mercator (projection) 

UV Ultraviolet 

VAX Virtual Address Extension 
WORM Write Once Read Many (times) 

Symbols 

a Formulation coefficient 
A(k) Absorptivity. 

b Formulation coefficient 
B Excess target radiance 

C Chlorophyll a pigment, or just pigment concentra- 
tion. 

[C + P] Pigment concentration defined as mg chlorophyll a 
plus phaeopigments m -3 . 

C13 Pigment concentration derived using CZCS bands 1 
and 3. 

C23 Pigment concentration derived using CZCS bands 2 
and 3. 

di Distance from the zth observation point to the point 
of interest. 

dj Distance from the jth observation point to the point 
of interest. 

Pfcm Percentage of energy removed from a wavelength 
band. 

Fo(\) Incident solar irradiance. 

G( A) P a (AO/Pa (670)== (670/A) 7 r 2r (670)/T 2r (Ai) 

i Interval index. 

/o Incident radiant intensity. 

/ 1 Radiant intensity after traversing through an ab- 
sorbing medium. 

I2 Reflected radiant energy received by the satellite 
sensor. 

j Interval index. 

k Wavenumber of light (1/A). 
k\ Beginning wavenumber 
/c 2 Ending wavenumber 

L(A) Radiance. 

L a (A) Aerosol radiance. 

L c ( A) Cloud radiance threshold 
L r (A) Rayleigh radiance. 

L r0 (A) Rayleigh radiance at standard atmospheric pressure, 
Po- 

L s ( A) Subsurface water radiance. 

L t (A) Total radiance at the sensor. 

L U (X) Upwelled spectral radiance. 

Lw{ A) Water-leaving radiance. 
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m Index of refraction. 

n(A) An exponent conceptually similar to the Angstrom 
exponent. 

N Total number of points or pixels. 

P Phaeopigment concentration. 

Po Standard atmospheric pressure (1,013.25 mb). 

P a Phaeopigment concentration. 

q Water transmittance factor. 

R Reflectance. 

R a Aerosol reflectance. 

Ra Ra/(qT 2r )- 
R r Rayleigh reflectance. 

R s Subsurface reflectance. 

Rt Total reflectance at the sensor. 

R t (R t - R r )/(qT 2r ). 

S( A) L a (A)/L a (670). 

t\ First observation time. 
t 2 Second observation time. 

t(k) Spectral transmission as a function of wavenumber. 
t( A) Diffuse transmittance of the atmosphere. 

T 2r Two-way diffuse transmittance for Rayleigh atten- 
uation. 

T oz Diffuse transmittance of ozone. 

Vi(tj) The ith spatial location at observation time tj. 

W Wind speed. 

x Abscissa or longitudinal coordinate. 

y Ordinate or meridional coordinate. 

a Percent albedo, tilt angle, or formulation coefficient 
(intercept). 

f3 Formulation coefficient (slope) 

7 Angstrom exponent. 

Ak Equivalent bandwidth. 

A P Pressure deviation from standard pressure, Po. 

6 Satellite zenith angle. 

0o Solar zenith angle. 

On The zenith angle of the vector normal to the surface 
vector for which glint will be observed. 

A Wavelength of light. 

p Mean value or cosine of the satellite zenith angle. 
fio Cosine of the solar zenith angle. 

Vj The jth temporal weighting factor. 

p Fresnel reflectivity. 

cr Standard deviation. 

a 2 The mean square surface slope distribution. 

T 0 x Oxygen optical thickness at 750 nm. 
r r Rayleigh optical thickness. 

r r o Rayleigh optical thickness at standard atmospheric 
pressure, Po- 

Satellite azimuth angle. 

4>o Solar azimuth angle. 

u; The surface reflection angle. 
cji Spatial weighting factor. 
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COLOR PLATES 
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Plate 1 . Pigment concentration images derived for scene 1 using the Gordon 2-band, Clark 3-band, 
Smith-Wilson, and European methods. 
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Plate 2. The Angstrom exponents derived from the European method for scene 1. 
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Plate 3. Pigment concentration images derived for scene 2 using the Gordon 2-band and European 
methods. 
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Plate 4. The Angstrom exponent derived from the European method for scene 2. 
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Plate 5. Composite of three pigment concentration images using three different ozone concentrations 
and the frequency distributions of log [pigment concentration]. The color ranges are 0-0.10 (purple hues), 
0.10-0.45 (blue hues), 0.45-1.0 (green hues), 1.0-1. 6 (yellow hues), 1.6-3.0 (orange hues), 3.0-4. 5 (red 
hues), greater than 4.5 (brown hues). Mean pigment concentrations for the entire image and the ozone 
concentration are given at the top of each panel. 
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Plate 6. Composite of three pigment concentration images using three different surface pressures and 
the frequency distributions of log[pigment concentration]. The color ranges are 0-0.10 (purple hues), 
0.10-0.45 (blue hues), 0.45-1.0 (green hues), 1.0-1. 6 (yellow hues), 1. 6-3.0 (orange hues), 3. 0-4. 5 (red 
hues), greater than 4.5 (brown hues). Mean pigment concentrations for the entire image and the surface 
pressure are given at the bottom of each panel. The fourth panel shows the frequency distributions of 
log [pigment] for the three analyses. 
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Plate 7. Pigment concentration image derived using center of scene ozone and standard surface pressure 
values. 
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Plate 8. Pigment concentration image derived using a pixel-by-pixel pressure correction and the center 
of scene ozone value. 
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Plate 9. Pigment concentration difference image, Plate 8 minus Plate 7. 
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Plate 10. Pigment concentration image derived using a pixel- by- pixel ozone correction and a standard 
pressure value. 
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Plate 11. Pigment concentration difference image, Plate 10 minus Plate 7. 


ORIGINAL PAGE 

COLOR PHOTOGRAPH 





C. McClain, J. Comiso, R. Fraser, J. Firestone, B. Schieber, E. Yeh, K. Arrigo, and C. Sullivan 



Plate 12. Pigment concentration image derived using both pixel- by- pixel pressure and ozone correc- 
tions. 
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Plate 13. Pigment concentration difference image, Plate 12 minus Plate 7. 
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Plate 14. Unmapped (upper two panels) and mapped (lower two panels) pigment concentration images. 
The color ranges are 0-0.10 (purple hues), 0.10-0.45 (blue hues), 0.45-1.0 (green hues), 1.0-1. 6 (yellow 
hues), 1. 6-3.0 (orange hues), 3.0-4. 5 (red hues), greater than 4.5 (brown hues). 
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Plate 15. Monthly average total ozone images, based on data from NIMBUS-7 TOMS collected be- 
tween November 1978 and January 1992. The images were generated with SDSIMAGE and SDSNANIMATE, 
procedures written using IDL. 
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Plate 16. Monthly average surface wind speed images, based on data from COADS between 1946 and 
1990. The images were generated with SDSIMAGE and SDSANIHATE, procedures written using IDL. 
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Plate 17. Similar to Plate 16, except illustrates monthly average surface pressure images. 
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Plate 18. Similar to Plate 16, except illustrates monthly average surface relative humidity images. 
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PLATE 19. Pigment concentration images using no ringing mask, the SEAPAK mask (DISTANCE— 20), 
the Mueller mask, and the Miami DSP mask generation methods. 
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PLATE 20. Unmasked (upper left panel) and sun glint flagged (upper right and lower two panels) pigment 
concentration images. Color ranges are 0-0.10 (purple hues), 0.10-0.45 (blue hues), 0.45-1.0 (green hues). 
Wind speeds to generate the glint flags are 1 (pink), 3 (red), 5 (green), 7 (yellow), 9 (orange), 11 (cyan), 
and 13 ms" 1 (tan). 
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